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I. INTRCDOCTION 

Simulations of observed tropical climate by general circulation 
models (QCMs) have steadily improved over the last three decades as 
refinements have been made in the treatment of interactive physical 
processes and in the numerical techniques employed to solve the 
governing equations. Not only do we know that the tropics is a major 
energy source which drives the atmospheric motions, but we new need to 
examine the intricate interaction between middle latitudes and 
tropics. This investigation may have an inportant bearing on our 
understanding of long-term fluctuations of the global climate. 

The large-scale motions in the tropics, which are quite different 
from those of midlatitude systems due to baroclinicity of the 
atmosphere, are driven primarily by latent heat released from strong 
cumulus convections. Because of the special nature of this driving 
force as well as the smallness of the Coriolis parameter, the large- 
scale motion system has certain distinctive structural features such 
as monsoon circulation associated with the tropical region of Asia and 
Africa, and east-west Walker circulation spanning the breadth of the 
Pacific Ocean. 

Based upon a shallow water model, Matsuno (1966) investigated the 
nature of large-scale waves in the equatorial area. He found that 
two fundamental large-scale wave modes could exist in the equatorial 
atmosphere: westward moving mixed Rossby-gravity waves and eastward 

moving Kelvin waves. Later, observational studies by Yanai and 
Maruyama (1966) and Wallace and Kousky (1968) confirmed the existence 
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of the equatorial planetary waves. These discoveries revived interest 
in wave disturbances in the tropics and have prompted subsequent 
observational and theoretical studies on this topic [see review 
articles by Wallace (1973) and Holton (1975) ] . At present, three 
ideas appear true. First, both the Kelvin waves and the mixed Rossby- 
gravity waves are excited in the upper troposphere. Second, they are 
the predominant disturbances of the equatorial stratosphere. And 
third, through their vertical energy and momentum transport, the waves 
play a crucial role in the general circulation of the stratosphere. 

In addition to their role in the stratospheric quasi-biennial 
oscillation, equatorial planetary waves also are believed to be 
responsible for the eastward propagating 30-50 day oscillations. 

Madden and Julian (1971; 1972) first described the 30-50 day 
oscillation as global-scale eastward propagating zonal circulation 
cell along the equator. Due to the propagating feature and time scale 
of the oscillation, Chang (1977) suggested that the eastward 
propagation is the manifestation of convectively driven equator ically 
trapped Kelvin waves. Following the same line, Lau and Peng (1987) 
recently proposed that the observed eastward propagation of 
intraseasonal oscillation in the tropical atmosphere arises as an 
intrinsic mode of oscillation resulting from an interaction of 
convection and dynamics via the so-called "mobile" wave-CISK 
mechanism. 'Through this mechanism, the heat source feeds on the east- 
west asymmetry of forced equatorial waves. As a result, Kelvin waves 
are selectively amplified and, in turn, cause the heat source to 
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propagate eastward. 

therefore, it is essential to simulate the equatorial planetary 
wave in a proper manner with a GCM before attempting to use it for 
climatic studies. Equatorial planetary waves have been simulated by 
both the GCM of the Geophysical Fluid Dynamics Laboratory (GFDL) 
(Hayashi, 1974) and by the first generation GCM of the National Center 
for Atmospheric Research (NCAR) (Kasahara et al. , 1973) . Since the 
NCAR Community Climate Model (CCM1) is the latest version of the 
spectral GCM at NCAR, we will examine its ability to successfully 
simulate the equatorial planetary waves. We shall also examine how 
these waves are maintained. 

There have been several theoretical studies to explain the 
generation and maintenance of the equator ial planetary waves (see 
Chapter III for details) . In view of the controversies among the 
proposed linear theories, we feel that numerical experiments with a 
GCM may be used to clarify same issues on the generation of the 
tropical disturbances. Three experiments with the NCAR (XML under 
seasonal variation conditions were analyzed in this study. The first 
is a control experiment run with complete physics. The second and the 
third experiments run without the effects of topography and clouds, 
respectively (see Chapter III for details) . 

The outline of this study is as follows: Chapter II reviews the 

past studies on the equatorial planetary waves. Chapter III describes 
the basic GCM used in this study and the three experiments. The data 
and computational schemes used in this study are discussed in Chapter 
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IV. The results are presented in Chapter V. Chapter VI presents the 
concluding remarks and future studies. 
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II. HISTORICAL DEVELOPMENT 

There has been significant progress in observational and 
theoretical studies of equatorial planetary waves since the discovery 
of mixed Rossby-gravity waves and Kelvin waves in the equatorial 
stratosphere. In this chapter, we shall review the past studies on 
equatorial planetary waves. The review consists of three sections: 
Theoretical Work, Observational Studies, and Model Simulations. 

A. Theoretical Work 

The rotation of the earth dominates the way the atmosphere and 
ocean respond to imposed changes. This dynamic effect is caused by 
the Coriolis acceleration, which is equal to the product of the 
Coriolis parameter and the horizontal velocity. So the adjustment 
processes are somewhat special when the Coriolis acceleration 
vanishes. 

Oceanographers found that an ocean channel enclosed in an 
equatorial zone can have an oscillatory mode whose period is much 
larger than one day. Stem (1963) first suggested this idea by 
considering inertial gravity waves in a rotating, homogeneous fluid 
confined between two concentric spheres of nearly equal radii. He 
shewed that there exist standing oscillations of frequency a( xib/R)^^, 
where o is the rate of rotation of the system, h is the thickness of 
the fluid, R is the mean radius of the spheres, and n is a positive 
integer. His analysis assumed that b/R « 1, and he showed that the 
wave amplitude decreases to zero outside a band of order (h/R)^^ on 
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either side of the equator. 

Bretherton (1964) gave a sinple explanation of this trapping 
phenomenon in terms of wave reflection. He suggested that outside a 
certain band, the Coriolis parameter is greater than the frequency of 
waves inside that band. For the waves in the band, both gravitational 
and Coriolis forces tend to restore a particle to its original 
latitude. Thus, there is a minimum frequency for the oscillation at a 
given latitude, namely, that of purely horizontal motions under the 
action of the Coriolis force. A wave travelling northward or 
southward must conserve this frequency, and cannot penetrate into a 
region where the Coriolis parameter exceeds it. So the wave is 
reflected and a standing oscillation is possible. He discussed 
baroclinic internal waves in the equatorial ocean using a two- layer 
model, and he showed that one such oscillation would have a period of 
about 10 days and would be trapped between 6N and 6S. 

In order to study the atmospheric motions in the equatorial area, 
Matsuno (1966) derived a dispersion relation for the wave solution of 
the shallow water equations on an equatorial /3-plane and gave the 
horizontal depictions of the wave disturbances, which are a set of 
equatorially trapped wave modes. Later, the theory of equatorial 
waves was extended to an isothermal atmosphere by Lindzen (1967) and 
Holton and Lindzen (1968) . 

The resulting family of the wave mode in the shallow water model 
contains gravity and Rossby waves. The gravity waves propagate both 
eatward and westward, while Rossby waves only propagate westward with 
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slow velocity. These two types of waves are distinguished from each 
other by the difference of their frequencies. Namely, frequencies of 
the gravity waves are much larger than those of the Rossby waves. 

Figure la shews the horizontal pattern of the westward moving 
wave associated with the lowest mode in the wave family. This wave, 
with the smallest north-south extent, has antisymmetric structures in 
the temperature and zonal wind, and a symmetric structure in 
meridional wind with respect of the equator. Examining this figure, 
we note that this wave has somewhat mixed characteristics of the 
Rossby and gravity waves. The relationship between the pressure and 
the velocity fields is approximately geostrophic in the higher 
latitudes. But near the equator, ageostrophic wind components 
dominate. The configurations of wind and pressure fields near the 
equator resemble those of the gravity wave. The dispersion 
relationship of this wave is similar to that of the Rossby wave when 
the zonal wavelength is smaller than the meridional extent. For the 
long-wave part, however, the frequency becomes very long and 
approaches that of the gravity wave of the same wavelengh. For this 
reason, the wave of this type is called a mixed Rossby-gravity wave. 

By assuming zero meridional velocity, the shallow water equation 
has a special solution, which corresponds to an eastward moving wave 
with a phase velocity of prue gravity wave. Figure lb shows the 
horizontal structure of this wave. It has symmetric structures across 
the equator in pressure and zonal wind fields and has no meridional 
velocity component. Figure lb also shows that the wave behaves like a 
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pure gravity wave in zonal direction. But in the north-south 
direction the geostrophic relation holds between zonal velocity and 
meridional pressure gradient, Hiis type of wave is called a Kelvin 
wave which propagates along a coastal boundary and has no velocity 
component normal to the boundary. In the atmospheric case, the 
equator plays the same role as the coast line. 

Figure lc shows the horizontal distribution of the pressure and 
velocity fields of the lowest mode Bossy wave. Hie wave is 
characterized by the geostrophic balance between pressure and velocity 
fields. It is interesting that strong zonal velocity is found along 
the equator. We expect this strong velocity from the approximate 
geostrophic balance between the pressure gradient and wind. 

A very important property of the equatorial disturbances is 
trapping of the lcw-frequency waves. In other words, the disturbances 
are confined to the vicinity of the equator. This trapping effect is 
entirely due to the variation of the Coriolis parameter. In the case 
of the inertial gravity wave, the trapping phenomenon can be 
interpreted from the viewpoint of reflection of the primary wave as 
explained by Bretherton (1964) . Hie wave generated near the equator 
will be refracted and reflected toward the equator. 

On the other hand, in the case of equatorial Kelvin wave, which 
propagates eastward with the velocity of a long gravity wave, the 
deviation of the surface elevation is just in geostrophic balance with 
the motion associated with the wave. Also, the amplitude of the wave 
decreases exponentially away from the equator. Hiis changing of the 
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Coriolis parameter's sign near the equator confines Kelvin waves only 
in the vicinity of the equator. Here, the equator plays the role of a 
duct in the propagations of these waves. 

Lindzen (1967) studied the vertical structures of equatorial 
waves in an isothermal atmosphere. He found that waves with periods 
of several days, which generally decay with height in the middle 
latitude, may become untrapped and propagate vertically near the 
equator due to the decreasing Coriolis frequency. 

Vertically propagating wave disturbances in the equatorial area 
influence the atmospheric general circulation by transporting energy 
and momentum upward. Holton and Lindzen (1968) discussed the 
characteristics of the vertically propagating Kelvin wave in the 
tropical atmosphere. They found that the Kelvin wave with eastward 
tilting phase line is likely to be excited in the upper troposphere 
and transports westerly momentum upward. Therefore, Lindzen and 
Holton (1968) proposed a mechanism whereby interactions between the 
wave disturbances and the mean zonal flow could provide the momentum 
source for the quasi— biennial oscillation. On the other hand, the 
mixed Rossby-gravity wave, which tilts westward with height and 
transports westerly momentum upward, may be partly responsible for the 
easterly acceleration of the quasi-biennial oscillation of the 
equatorial stratospheric zoned, wind. 

The vertically propagating waves are continuously carrying energy 
upward from the lower atmosphere to an energy sink at a higher level. 
The energy sink may be the mean zonal flow. Or it may be eddy 
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viscosity or radiative dumping, if the waves propagate to very high 
levels without being absorbed by the mean zonal flow. Either case 
implies a source of wave energy in the lower atmosphere. 

To shew the steady circulation in the shallow water model, 

Matsuno (1966) worked out an example that included mass source and 
sink imposed altematingly along the equator. Hie remarkable feature 
of the resultant motion is the strong zonal flow confined in the 
vicinity of the equator. This flew ridge and trough are located along 
the equator, which divides the pressure cell into two petals. 
Mathematically, this result can be explained by the fact that Rossby 
wave solution predominates while gravity waves are suppresed. 

B. Observational Studies 

Theoretical contributions by Matsuno (1966) and Lindzen (1967) 
have provided invaluable guidance in the planning of synoptic 
investigations and the interpretation of results. The existence of 
the mixed Rossby-gr avity wave was confirmed by Yanai and Maruyama 
(1966) from analyses of wind data in the lower stratosphere over the 
Marshall Island in the equatorial Pacific. They noticed regular, 
wavelike fluctuations in the meridional wind component at 
stratospheric level. These disturbances were observed to propagate 
westward and downward with a horizontal wavelength of about 10,000 km, 
a vertical wavelength of about 6 km, and a period of 4-5 days. 

Subsequent studies (Maruyama and Yanai, 1967; Maruyama, 1967) 
revealed the existence of associated fluctuations in the zonal wind 
component and in temperature. The observed phase relationships rotate 
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clockwise with time as waves pass, and the maximum temperature occurs 
at the time of the maximum southerly wind. Hie reverse is true for a 
point in the Southern Hemisphere. Hie waves have antisymmetry about 
the equator and appear to be confined within 12° of the equator. Hie 
horizontal structure and phase velocity of the observed disturbance 
showed close agreement with those of the mixed Rossby-gravity wave 
discussed by Matsuno (1966) . 

Lindzen and Matsuno (1968) further studied the mixed Rossby- 
gravity wave, obvserved by Yanai and Maruyama (1966) , in terms of 
forced oscillations. Assuming the wave was excited in the lower part 
of the atmosphere, they were able to explain the importance of the 
wave's tilts and the wave's contributions to the mean zonal momentum. 

Yanai et al. (1968) used the cross-spectrum analysis method to 
describe the more prominent fluctuations in the meridional wind 
component at tropospheric levels. Hie researchers found a wave mode 
at levels above 250 mb with a period of 4-5 days. Hie mode displays 
the same characteristics as the stratospheric mixed Rossby-gravity 
wave observed by Yanai and Maruyama (1966) . Yanai et al. also 
established that the mixed Rossby-gravity waves exist in both upper 
troposphere and lower stratosphere. 

Hie analysis of the large-scale wave in the lower stratosphere 
suggested the importance of inquiries into the upper troposphere, 
where the vertically propagating waves might be excited. Maruyama 
(1968b) and Yanai and Hayashi (1969) calculated vertical and 
horizontal transports of energy in the upper troposphere and the lower 
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stratosphere. They obtained the following results: 

1. Hie vertical axis of the large-scale equatorial waves that 
move west relative to the basic wind have a westward tilt. 
That tilt indicates that the wave energy is transported 
upward, if the wind shear is not too strong. 

2. The large-scale waves in the equatorial upper troposphere 
provide the lower stratosphere with wave energy. The 
vertical convergence of energy occurs at levels from 15-20 
km. 

3. In the lower stratosphere, horizontal convergence of wave 
energy occurs near the equator. 

Soon after the discovery of 4-5 day stratospheric waves, Wallace 
and Kousky (1968) provided observational evidence for the existence of 
the disturbances in the zonal wind component with eastward and 
downward propagating phase line in the equatorial lower stratosphere. 

A distinguishing characteristic of the wave is the absence of 
fluctuations of the meridional wind. The zonal wind and pressure 
field were found to be in phase and oscillate with a period of 10-15 
days and zonal wavenumber 1-2 with vertical wavelength of 5-10 tan. 

This oscillation was identified as an atmospheric Kelvin wave 
discussed theoretically by Matsuno (1966) , Lindzen and Holton (1968) , 
and Holton and Lindzen (1968) . 

The observed temperature field of a Kelvin wave has a distinct 
fluctuation. Warmest temperatures precede the maximum westerly by 1/4 
cycle, with both oscillations propagating downward. The wave 
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amplitude is largest at the equator and decays to about half the 
maximum value at 10° latitude. 

This observational evidence suggests some physical 
interpretations as follows: 

1. The absence of meridional wind fluctuations requires that the 
zonal wind component must be in geostrophic equilibrium with 
meridional pressure gradient. It follows that the levels of 
maximum zonal wind and maximum pressure should coincide. 

2. Hydrostatic equilibrium requires that the warmest temperature 
be 1/4 wavelength below the level of the maximum pressure. 

3. Hie zonal momentum equation requires that maximum westerlies 
follow 1/4 cycle behind the maximum east-west pressure 
gradient. 

4. In the absence of the meridional wind component, continuity 
requires that the zonal wind and the vertical motion 
oscillations be in phase. 

Hie upper stratospheric wave with perod of 5-10 days found by 
Hirota (1978, 1979) is also likely to be a Kelvin wave and is 
associated with a larger vertical wavelength (20-25 km) than the lower 
stratospheric Kelvin wave. This upper stratospheric Kelvin wave may 
play an important role in the semiannual oscillation in the tropical 
middle atmosphere through a possible upward transport of eastward 
momentum as suggested theoretically by Dunkerton (1979; 1982) . 

At mesospheric levels, there is also evidence of faster eastward 
propagating disturbances. Salby et al. (1984) reported fluctuations. 
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identified with Kelvin mode associated with wavenumbers 1 and 2, in 
satellite observation between 16 and 70 km. The wavenumber 1 
component had an eastward moving period of 35-40 days, and an eastward 
tilt of 40 km vertical wavelength, while the wavenumber 2 component 
had an eastward moving period of 3. 8-4 .3 days, and an eastward tilt of 
about 20 km vertical wavelength. 

Recently, Zangvil and Yanai (1980) , and Yanai and Lu (1983) 
presented a wavenumber-frequency analysis of the large-scale transient 
wave disturbances. The disturbances were at the 200 mb level in the 
latitude belt from 20S to 40N. The analysis was based on wind data 
collected during two northern simmers, 1967 and 1972. Decomposing the 
wind data into symmetric and antisymmetric components with respect to 
the equator (Yanai and Murakami, 1970) and using the wavenumber- 
frequency analysis technique developed by Kao (1968) and Hayashi 
(1971; 1973) , they were able to separate the equatorial disturbances 
into eastward and westward moving modes. Three types of equatorial 
waves were detected in the wind data of 1967 ; the Kelvin wave 
associated with zonal wavenumber m=l and 2 and with a period of 7 
days, and itf=1 with a period longer than 20 days; a mixed Rossby- 
gravity wave associated with zonal wavenumber 4 and a period of 5 
days; and an equatorially trapped Rossby wave with m=2 and a period 
near 10 days. On the other hand, in 1972, the mixed Rossby-gravity 
wave was absent and only weak signals of Kelvin waves with m=l and 2 
and the equatorally trapped Rossby wave with m=2 were found. 

The researchers also showed that all the equatorially trapped 
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westward moving waves are associated either with equatorward energy 
flux at both 3 ON and 15S or with strong equatorward energy flux at 3 ON 
and weak poleward energy flux at 15S. They also showed that clear 
correspondence of the major peaks of meridional convergence of wave 
energy flux in the tropical belt to the equatorially trapped westward 
moving waves. Their findings strongly suggests that the westward 
moving waves may be excited by presure work due to mid-latitude 
disturbances, i.e. , lateral forcing. However, there was no evidence 
to suggest the association of the wave energy flux convergence with 
Kelvin waves. 

Zangvil and Yanai (1980) analyzed cloud activities in the 
latitude belt 20S and 40N for the period June-August 1967 and found 
cloud activity was accompanied by the long-period eastward moving 
waves, whose characteristics were quite consistent with those of 
Kelvin waves. 

In order to compare our results with the observational studies, 
we summarize in Table 1 the properties of Kelvin and mixed Rossby- 
gravity waves that have been identified observationally. 

C. Model Simulations 

Remarkable progress has been made in the linear theory and the 
analysis of the equatorial planetary waves. In order to have a better 
understanding of the tropical disturbances, the characteristics of the 
waves appearing in the general circulation models at both NCAR and 
GFDL have been analyzed. The features of the two models are 

| 

summarized in Table 2. 

i 
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Table 2. Summary of three GCMs uesd to simulate equatorial planetary waves 
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Kasahara et al. (1973) and Tsay (1974) examined the 
characteristics of vertically propagating equatorial planetary waves 
appearing in a January simulation of the NCAR 12-layer stratospheric 
global circulation model with the top of the atmosphere extending to 
36 km. Three types of dominant waves were identified by spectral 
peaks and phases in the model's lower stratosphere: the Kelvin waves 

with zonal wavenumber m=l and 2, and with eastward-moving periods of 
24 and 12 days, respectively; the mixed Rossby-gravity wave with zonal 
wavenumber m=2 and a westward-moving period of 4.6 days; and the 
Rossby wave type of free oscillation with zonal wavenumber m=l and a 
westward-moving period of five days. Both Kelvin and mixed Rossby- 
gravity waves have a vertical wavelength of about 8.5 km. The 
amplitude of the Kelvin wave is about the same order of magnitude as 
that observed in the real atmosphere. However, the amplitude of the 
mixed Rossby-gravity wave is one order of magnitude smaller than that 
in the real atmosphere. The reason for this may lie with the 
prevailing easterlies in the model lower stratosphere. 

Tsay (1974) showed that the Kelvin wave probably was determined 
by a certain filtering effect of the atmosphere since no peak in the 
tropospheric heating function existed at the wavenumber and 
frequencies of the Kelvin wave. Tsay also showed that the possible 
energy source of the Kelvin wave is the meridional energy convergence 
between 12 and 18 km over the equator. 

The characteristic of the simulated Kelvin waves agreed well with 
those of observed Kelvin waves. However, the simulated wavenumber 1 
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oscillation had a significantly longer period. The vertical 
resolution or the forcing characteristics of the model would be 
responsible for the difference. Also, observations do not indicate 
wavenumber 2 Kelvin wave with period one-half that of wavenumber 1 
oscillation. On the other hand, such a wave was found in Holton's 
(1973) linearized model with white-noise forcing. They theorized that 
the frequency was determined by the vertical wavelength that best fit 
the vertical scale of the forcing field. 

Following the same line, Hayashi (1974) analyzed the equatorial 
waves simulated by the GFDL 11-layer general circulation model (for 
details, see Table 2) . He included seasonal variation during July 
through October. In addition to easterly waves, three types of 
traveling waves were isolated from stationary waves and ultralong 
waves extending from middle latitudes: 

1. Kelvin waves with wavenumber 1-2 and a period of 15 days. 
These waves attain their maximum amplitudes at 38 mb with 
their phase lines tilting eastward and westward with height 
in the stratosphere and the troposphere, respectively. 

2. Mixed Rossby-gravity waves with wavenumbers 3-5 and a period 
of 3-5 days. These waves attain their maximum amplitudes at 
110-190 mb with their phase lines tilting westward and 
eastward in the stratosphere and the troposphere, 
respectively. 

3. Rossby-type waves with wavenumbers 1-5 and a period of 15 
days. These waves attain their maximum amplitudes at 110-190 
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mb with their phase lines tilting westward and eastward in 
the stratosphere and the troposphere, respectively. 

All of these traveling waves are maintained primarily by the 
conversion of available potential energy generated by condensational 
heating. Due to the coupling between large-scale disturbances and 
convection, presumably through moisture convergence as in the case of 
the CISK theories, a broad spectrum of waves are generated and 
maintained by the latent heat released by the rainfall. Hie spectra 
are not random noise, but instead, exhibit certain characteristic 
wavenumbers and phase velocities. 

Hay ash i et al. (1984) also analyzed the large-scale equatorial 
disturbances simulated by the 40-level GFEL general circulation model 
with annual mean condition. Even though the model had no seasonal 
cycles, the large-scale equatorial eastward moving disturbances were 
still successfully simulated and agreed with observations by Wallace 
and Kousky (1968) , Hirota (1978, 1979) , and Salby et al. (1984) . 
However, the model did not detect mixed Rossby-gravity waves. 

There are, in general, three theories to explain the generation 
and selection of large-scale equatorial disturbances. The first 
approach is the theory of resonance to midlatitude forcing (Mak, 

1969) . Mak found that westward moving equatorial waves exhibit 
a large resonant response at a period of four days with wavenumber 4. 
He also found that the waves derive their energy from midlatitude 
disturbances. His result is consistent with wave propagation 
theories. According to Charney (1969) , waves moving westward relative 
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to the basic wind can propagate into the tropics, provided the 
horizontal shear is not too strong and the vertical scale is not too 
short. 

The second approach was an instability (CISK) theory by Hay ash i 
(1970) . The theory assumes that convective heating is proportional to 
surface convergence due to large-scale waves. The researchers found 
that the instability allows the existence of mixed Rossby-gravity 
waves and Kelvin waves with realistic three-dimensional structure. 
However, waves preferred scale did not agree with their observed. 

This difficulty can be avoided, however, by assuming that cumulus 
heating becomes inefficient for high-f requency waves due to the 
constraint of moisture availability (Kuo, 1975) . Moreover, the CISK 
theories also predict sharp spectral peaks in the tropical rainfall 
which are not clearly detected in general circulation models 
(Tsay, 1974; Hayashi, 1974) . 

The third approach is to examine the response of the atmosphere 
to random thermal forcing in the presence of a vertical wind shear, 
since the atmospheric heating field does not have spectral peaks at 
the wavenumbers and frequencies of the observed Kelvin waves and mixed 
Rossby-gravity waves. Holton (1972) demonstrated that both eastward 
and westward moving waves were excited in the troposphere by a 
localized tropospheric heat source pulsating with an imposed 
periodicity. However, depending on the vertical profile of wind, 
eastward moving Kelvin waves or westward moving mixed Rossby-gravity 
waves propagate into the stratosphere. Holton (1973) further pointed 
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out that a 10-30 day spectral peak in tropospheric thermal forcing is 
not required to account for the observed spectral distribution of 
stratospheric Kelvin waves, since the atmosphere acts as a bandpass 
filter to select these waves. This point was confirmed by the 
spectral analysis of the general circulation model (Tsay, 1974; 

Hayashi, 1974). 

In view of the controversies among these linear theories, Hayashi 
and Golder (1978) made an attempt to clarify, to same extent, the 
causes of the generation of an equatorial planetary wave. He conducted 
control experiments with a GFDL 13-level general circulation model (for 
details, see Table 2) . By successively eliminating the effects of 
topography, midlatitude disturbances, and condensational heat from the 
model, he showed the following results: 

1. The amplitude and the characteristic scale and period of both 
Kelvin and mixed Rossby-gravity waves are hardly affected by 
the land-sea contrast or the zonal variation of sea surface 
temperature. 

2. Kelvin waves are hardly affected by midlatitude disturbances. 
But the waves diminish drastically after condensational heat 
is also eliminated. 

3. Mixed Rossby-gravity waves are significantly reduced in 
amplitude by the elimination of midlatitude disturbances. 

They diminish completely after latent heat is also 
eliminated. 

4. The total kinetic energy of all equatorial transient 
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disturbances is generated primarily as the result of latent 
heat release in the tropics. That of subtropical transient 
disturbances is reduced significantly in the absence of 
midlatitude disturbances. 

It is of interest to compare Hayashi's results with the existing 
theories mentioned above. His results are consistent with these 
theories in that the topographical effect is not an essential factor 
for the charateristic scale and period of both Kelvin and mixed 
Rossby-gravity waves. Ihe results are similar to those of both the 
CIKS and random thermal forcing theories in that both Kelvin and mixed 
Rossby-gravity waves appear in the stratosphere as a result of latent 
heat release in the troposphere. On the other hand? the present 
results are also similar to those of lateral forcing theories in that 
tniv<vi Rossby-gravity waves can be significantly enhanced by 
midlatitude disturbances. However, these waves are not generated by 
midlatitude disturbances in the absence of condensational heating. 

Although the above theories have their own drawbacks, they should 
not be entirely rejected since the nonlinear effect is neglected. 
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III. MODEL AND EXPERIMENTS 
A. The Model 

The GCM used in this study is the latest version of the NCAR 
Community Climate Model version 1 (CCM1) . This version has evolved 
over the last few years from the earlier CCMOB. Substantial changes 
have been incorporated into the radiation scheme in both the physical 
parameterizations and the actual structure of the code. 

This model was chosen because it is designed for global forecast 
studies (both medium- and long-range) as well as climate simulations, 
and because it is widely uesd for both purposes by university as well 
as NCAR scientists. 

1_. The governing equations 

The governing equations are similar to those used by Bourke 
(1974) and Hoskins and Siimons (1975) . These equations adopt the 
vertical coordinate a = P/P s (p s is the surface pressure) proposed by 
Phillips (1957) . In the a -coordinate, the ground is a coordinate 
surface so that the effect of mountains can easily be incorporated 
without introducing problems associated with uncentered horizontal 
differencing. The governing equations include the following: 

• the vorticity (n) and divergence (6) equations; 

• the thermodynamic equation; and 

• the moisture (q) forecast and continuity equations. The 
continuity equation is not used directly but, when integrated 
in the vertical, gives equations for the surface-pressure 
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tendency and for b . 

Included in the model are the following parameterized physical 
processes: convection and condensation (Manabe et al. , 1965); fluxes 
of sensible and latent heat following a bulk aerodynamic formulation 
with the assumption of no heat storage over land; and interactions with 
subgrid-scale motions through horizontal and vertical diffusion. 
Moisture in excess of 100% relative humidity is precipitated without 
evaporation of the condensation in intervening layers. The algorithm 
is essentially that used in CCMO and represented in Ramanathan et al. 
(1983) . 

2. The numerical algorithm 

The model uses a semi-finite difference schane in space and a 
central and semi- implicit scheme in time. The vertical and temporal 
aspects of the model are represented by finite-difference 
approximations , while the horizontal aspects are treated by the 
spectral-transform method. In general, physical parameterizations and 
nonlinear operations are carried out in grid-point space. Horizontal 
derivatives and linear operations are performed in spectral space. 

The horizontal representation of an arbitrary variable q is expanded 
in a truncated series of surface spherical harmonics: 

M N(m) 

Q(A, u) * z i Q m p m (y) e lmX t 3 * 1 ) 

m*-M n= |M| n n 

where 

v = sin<*>. 
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0 = latitude, 

A = longitude, 

M = highest Fourier wavenumber included in the east-west 

representation, and 

N(m) = highest degree of the associated Legender functions 
included in the north-south representation (can be a 
function of the Fourier wavenumber (m) ) . 

In this study, the surface spherical harmonics are truncated 
rhamboidally at wavenumber 15. Ihe associated Gaussian grid has 48 
points in longitude with a spacing of 7.5 degrees and 40 points in 
latitude with an approximate spacing of 4.5 degrees. Hie algebra 
involved in the spectral transformations has been presented in several 
publication (Daley et al. , 1976; Bourke et al. , 1977; Machenbauer , 
1979) . 

Finite differences are used in the vertical. The model levels, 
at which a assumes the values of 0.991, 0.926. 0.811, 0.644, 0.500, 
0.336, 0.189, 0.074, and 0.009, are identical to those chosen by 
McAvaney et al. (1978) and to those used in many GFDL model 
simulations. Cn the vertical discrete grid, all prognostic variables 
and geopotentials are carried at full-index ^-levels, while the 
vertical velocity a is carried by half- index a -levels. The vertical 
finite differences are a combination of two differences: 1) the 

difference adopted by the Australian Numerical Meteorological Research 
Centre (ANMRC) model (Bourke, 1974; Bourke et al. , 1977) , and thus by 
CCM0 and 2) those differences used in the first European Center for 
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Medium Range Weather Forecast (EC3WF) model (Bur ridge and Haseler, 
1977) . The full-index ^-values of grid levels are specified as input 
to the model. The half- index levels are given as the average of the 
adjacent full-index levels with the end values set to 0 and 1. 

The time differences are basically centered with the terms 
responsible for fast-moving gravity waves treated semi- implicitly 
(Hoskins and Simmons, 1975) . The vertical and horizontal diffusion 
terms are treated implicitly using a time-splitting procedure — 
vertical before the advection processes and horizontal after. After 
the time step is completed, the temperature and mixing ratios are 
convectively adjusted. And finally, a time filter is applied to the 
prognostic variables lnp g , t, 7, 5 , and q. The time filter was 
originally designed by Robert (1966) and later studied by Asselin 
(1972) . 

3. The ini tied , and boundary conditions 

a. The initial conditions The basic initial data for the 
model consists of values of u, v, T, q, p g and 0 g on the Gaussian grid 
at time t=0. Fran these, U(=u cos#) , V(=v cos#), T' , and lnp_ are 
computed on the grid. Hie Fourier coefficients of these variables U 111 , 
V* T ,ra , lnp_ m are determined via an FFT subroutine. The 
spherical harmonic coefficients T' n m , q n m , lnP sn m , and 0 sn m are 
determined by Gaussian quadrature. Hie relative vorticity f and 
divergence 5 spherical harmonic coefficient are determined directly 
from the Fourier coefficients U m and V m using the relations 
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1 9V . 1 3U 
a(1-y 2 ) 8X a 3w 


(3.2) 


i au l 3V 

" a(l-y 2 ) 3 ^ ‘ 


(3.3) 


The first-time step of the model is forward semi-inplicit rather 
than centered semi- implicit. Therefore* only variables at t=0 are 
needed. 


b. The boundary conditions The required boundary data are 
surface temperature T g a t each ocean point, the surface geopotential 
at each point, the soil moisture at each ocean point, and the snow 
amount at nonocean points. However, these data are used only as the 
first guess in the T calculation the first time and are never used 
again. 

When the model is run with seasonal cycle, the sea-surface 
temperature (SST) , sea- ice locations, and snow cover are changed to 
reflect the seasonal changes. The SST and sea- ice distributions are 
specified from the monthly mean analysis of Alexander and Mobley 
(1976) . The monthly mean SST values are assigned the mid-month date 
and updated every 12 hours at each grid point using linear 
interpolation. 

The sea-ice and snow cover are updated at the beginning of every 
month and are held fixed during the month. The snow cover is a 
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function of latitude only. The radiation parameterization requires 
monthly mean ozone mixing ratios and the surface albedo. 

B. Model Simulations 

For the purpose of this study, three nunerical experiments with 
the NCAR CCM1 general circulation model were analyzed. The three 
experiments were the control experiment, the no-mountain experiment, 
and the no-cloud experiment. 

— 5?? control experiment 

Equatorial planetary waves have been well simulated by both the 
NCAR (Tsay, 1974) and GFDL (Hayashi, 1974) general circulation models, 
although the formulations of the GCM used by the two groups vary 
considerably. Since the CCM is now commonly used for climate studies 
for troposphere dynamics, the CCM1 is the latest version of the GCM 
with substantial changes incorporated in the radiation scheme. It is 
of interest to see to what extent the model can simulate the 
equatorial planetary waves. 

The control experiment, with various physics as described in the 
previous section, was initialized from real data in December 1, 1972 
and integrated for four years with seasonal forcing (i.e. , sea-surface 
temperature, sea-ice distribution and solar zenith angle) until the 
model reached its equilibrium state. Then, the experiment was run for 
an additional year. Only the last year's experiment instead of that 
of the first four years was used in order to avoid the effects of 
initial transients. 
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J2. Hie no-mountain experiment 

There have been numerous studies of the effects of large-scale 
orography on atmospheric flew. Hie major effects of large-scale 
orography which have been identified are: 

• The Siberian high appears concurrent with the appearance of 
the mountain and relates to the mountain's height. A 
stationary trough appears in the lees of the major mountain 
ranges (Manabe and Terpxstra, 1974; Mintz, 1965) . 

• Hie presence of mountains is instrumental in maintaining the 
South Asian low pressure system and helps to extend a monsoon 
climate farther north onto the Asian continent (Hahn and 
Manabe, 1975) . 

• Hie subtropical jet streams are weaker and the stationary 
waves have smaller amplitudes in the no-mountain case than do 
those in the mountain case (Blackmon et al. , 1987) . 

• Hie model without mountains produces a more zonally 
uniform distribution of tropical precipitation. 

Therefore, since there is no topographic effect to fix the 
stationary wave's nodes and antinodes, the no-mountain experiment 
weakens not only stationary waves, but also standing wave 
oscillations. Consequently, the interaction between midlatitudes and 
tropics, which is essentially contributed by the stationary ultralong 
waves, is also reduced. It will be of interest to examine whether 
reducing the interaction between midlatitude and tropics affects the 
characteristic scales and periods of the equatorial planetary waves. 
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In this experiment, all conditions remained the same as those in 
the control except for the elimination of topographical effects by 
replacing the model's surface with ocean points. The numerical time 
integration was then carried out with the same procedure as the 
control experiment. 

3_. The no-cloud experiment 

Cumulus convection and latent heat release play key roles in 
determining the structure of the temperature and moisture fields in 
the atmosphere. Heating and moistening, due to these processes, 
interact with radiation and large-scale dynamics in determining 
important aspects of the general circulation. Due to the monsoon type 
of the circulation, the energy is transported up f ran the warm sea- 
surface temperature area by means of cumulus convection and latent 
heat release in the tropics. This transport is confirmed by Donner et 
al. (1982) . Fran a GCM simulation, they showed that cumulus 
convection warms the upper troposphere and slightly cools the lower 
tropical troposphere. In general, the large-scale motions of the 
atmosphere organize the pattern of heating. According to the 
satellite observations of cloud cover compiled by Sadler (1970) , 
intense convection is associated particularly with the northern part 
of South ftnerica, equatorial Africa, and equatorial western Pacific. 
Also, the convection is associated with the intertropical convergence 
zone (ITCZ) and the South Pacific convergence zone (SPCZ) . One might 
expect that the heating due to the latent heat released by 
condensation is significant over these areas. Cumulus convection, as 
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feedback, influences the large-scale atmospheric motion both directly 
through heating and moistening and also indirectly by inducing changes 
in the meridional circulation. Chen and Baker (1986) calculated the 
global di aba tic heating during PGGE year and pointed out that the 
di aba tic heating in the tropics largely balanced by the vertical 
motion. 

Previous studies (Hayashi, 1970; Holton, 1972) suggested that the 
latent heat release, which is essential in both the CISK theories and 
thermal forcing theories, is the possible energy source of the 
equatorial planetary waves. A no-cloud experiment serves the 
purpose to examine the generation and maintenance of the equatorial 
planetary waves in the absence of cumulus convection in the 
atmosphere. 

The no-cloud experiment is identical to the control experiment, 
except that the cloud formation system in the model was turned off. 
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IV . DATA AND COMPDTATICNAL SCHEME 
A. Data 

In this study, we have analyzed three experiments of the NCAR 
CCM1 global general circulation model. The three experiments are 
control, no-mountain, and no-cloud experiments. For each experiment, 
we used two data samples during the periods December 1 to April 30 
(winter) and May 1 to September 27 (summer) , 1977, the fourth year of 
the model simulation. Daily simulations were stored for analysis. 

The data are gridded on an approximately 7 .5 ° x 4.5° mesh in the 
horizontal and seven levels in the vertical. Since each experiment 
was run continuously for five years, the data used in this study 
contain seasonal cycle. 

B. Computational Scheme 

This study involves several different statistical computational 
schemes. Among these schemes, space-time spectral analysis is 
estimated by the lag-correlation method. This section is devoted to a 
brief description of these schemes. 

1. Pre-filtering 

In order to isolate the tropical perturbation signals in which we 
are interested, we have designed a pre-filtering method which ronoves 
seasonal cycle and minimizes the signals with periods longer than 30 
days. 

Seasonal cycle is renoved by using the least squares fit method 
(Bendat and Piersol, 1971). Let {U }, n=l,2,...,N be the data values 
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sampled at an interval of h. Suppose it is desired to fit these data 
with a polynomial of degree k defined by 

K k 

U n = E ( nh ) n=l ,2, . . . ,N . (4.1) 

k-0 

In a least square fit, the set of coefficients { b^} are chosen so 
as to minimize the expression 

Q(b) = j, < U n ' * ", K - b k {„h) k 3 2 . (4.2) 

n-1 n=l k=0 

The desired coefficients can be obtained by taking partied 
derivatives of Equation (2) with respect to b and setting them equal 
to zero: 

an ' N K l 

IS-- E 2[U - l b k (r.hn[-(nh)*] . (4.3) 

dD l n=l n k=0 K 

This yields K + 1 equations of the form 

K N k+ , N 

E b k E (nh) k 1 = E U (nh)^ 
k=0 K n=l n=l n 

which now can be solved for the desired {b^} . 

Since we expect to remove the annual cycle mode in the data, we, 
therefore, choose k=2. In other words, we take second order 
polynomials to represent the seasonal cycle. Hie set of 
coefficients bgj b^ and b£ are determined by Equation (4.4) . 

In order to study the transient perturbances in the tropics, we 


1=0,1 ,2, . . . ,K 


(4.4) 
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nust separate these inodes from the lower frequency signals, which 
usually obscure the higher frequency signals because of larger energy 
content. This task is carried out by applying a second-order 
Butterworth bandpass filter fran the FDPACD package stored in the NCAR 
software library to the time series of the Fourier coefficients. By 
doing this, signals with periods smaller than 30 days are kept 
constant, while those with periods longer than 30 days are reduced. 


2. Lag-correlation and cross-power spectrum 

Following Bath (1974) , we assune a real time function f (t) having 
a fundamental period of 2” and satisfying Dirichlet's conditions can 
be represented by the following infinite Fourier series: 


f(t) - Z F(«) e" iwt 

ti)— —oo 


and F(u) is the Fourier coefficient given by 


F(u) 



f(t) e iujt 


dt . 


(4.5) 


(4.6) 


Now we consider the same conditions for another two time 
functions f-j (t) and ^2 ( t ) . The lag correlation function is defined as 
follows: 


+00 

Ci 2 ( r ) = / f-|(t) f 2 (t+r)dt 


(4.7) 


Cross-power spectrum is defined as 
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E ]2 (w) = F, (co) F 9 ( w ) 


(4.8) 


It can be proved immediately that the cross-correlation forms a 
Fourier pair with the cross-power spectrum: 


12 


< T > -h s 


+<» 


^1 2 ® 


iwr 


dio 


(4.9) 


■foo 


= C^Ct) e 1<i)r dr 


(4.10) 


In general, the cross-power spectrum E 12 (w)is complex: 


= + iQ-|2( a) ) 


(4.11) 


where k 12 (w) is the cospectrum and Q^ 2 (w) is the quadrature spectrum. 

3_. Space-time spectral analysis 

In order to interpret observed and simulated disturbances in 
terms of theoretical wave modes, it is helpful to analyze these 
disturbances into space-time Fourier components. However, due to the 
stochastic nature of atmospheric waves, these space-time Fourier 
components are randan variables and are not in themselves 
statistically and physically meaningful. For this reason, space-time 
spectral analyses (Kao, 1968; Hayashi, 1971; 1973) have been developed 
to obtain statistically meaningful results from a given space-time 
series which is cyclic in space and limited in time. 



38 


tet q(A,t) be a real, single-value function that is piecewise 
differentiable in a normalized domain, 0 < A,t £ 2v , where A and t 
stand for longitude and time, respectively. For convenience, a 
discrete representation is used for frequency as well as wavenumber. 
q(A,t) can be expanded into a space-time Fourier series as 

-|-oo *t*oo . z \ 

q(A,t) = l l Q(k,w) e l(kX+U)t) (4.12) 

k=-oo 

where Q(k,w) is the complex space-time transform determined by 

Q(k,u>) = -U/ 271 S™ q(A,t) e" l(kX+u,t) dAdt . (4.13) 

4iT 0 0 

Consider the same conditions for another scaler function q'(A,t) 
with a Fourier transform Q'(k,w). Thus, the contribution of 
and q’(A,t) integrated over a latitude circle and over a normalized 
time interval 2ir may then be expressed in terms of the sum of 
integrated over the frequency and wavenumber domain: 

i 2n 2 it +oo +oo 

— v S S q(A,t)q'(A,t) dAdt = l l Eqq * (k , +w) (4.14) 

4ir 0 0 k -o cj=o 

where Eqq' (k, to) is the cross-power spectrum of q(A,t) and q'(A / t) due 
to eddies of wavenumber k and frequency w. Hie positive and negative 
frequencies correspond to eastward and westward velocity, 
respectively. 

Hie cross-power spectrum between two sets of time series (q, q') 


| 
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is formulated as 

Eqq (k, +u>) = Kaj(C k ,C k ) + Kaj( » Sj^ ) 

± Q“(Ck,S k ) + Qo)(S k ,C k ) (4.15 

where K and Q denote cospectra, quadrature spectra, and C k and S k 
are the cosine and since space Fourier coefficients of q(A,t) . 

In particular, the power spectrum of the scalar quantity q(A,t) 
is given by 


Pq (k, jtu) = Eqq(k, +tu) 

= Pui(C k ) + Pa)(S k ) + 2Qto(C k ,S k ) (4.16) 

where Pw is the time power spectrum of a time series. 

A. Maintenance equation 

The maintenance of the equatorial disturbances will be 
investigated by using the energy transfer formulas derived by Hayashi 
( 1980 ). 

In spherical coordinates, the time change of wavenumber kinetic 
energy for large-scale motions is given by 

9K 

N n (k) - P n (., U ) 


sP n 4.v) sp„U,u) 

- [ — ■— — + — 1 

sy bp J 


( 4 . 17 ) 


where P n is the wavenumber cospectrum. Ihe kinetic energy per unit 
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mass, « , for the wavenumber n is defined by 


K n = I £ p n ( u >u) + P n (v.v)] , 


(4.18) 


and the nonlinear kinetic energy transfer spectrum N n (k) is defined by 

■ - tP„(u. $4 ♦ p n (v, |M1)] 


- [ p n(“- ff> + P n<*' 


+ t^ne [p^( u>uv ) . P n (v,uu)] 


- t p n< u - If 1 * P n^ v ’ lr )] 


(4.19) 


Here, 


ill = iU 

3y rcosex 


3( ) _ 3cos6( ) 

3y ~ rcos636 * 


The nonlinear energy transfer spectra N n (K) can be further 
partitioned into two parts as 


N (K) = <K„/K > + <K *K> . (4.20) 

n tti n on 

Here < K m 'K n > is the transfer of energy into wavenumber n by 
interaction among different wavenumbers excluding 0, while <K q - K p > is the 
transfer of energy into wavenumber n by interaction between the mean 
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flow and wavenumber n. 

By definition. <K -K > is given as 
J m n 


- TT-> + P „<*' TT-» 


3v'v' 


tane 


[P n (u, u V ) - P n (v, u'u' )] 


- C p n<“’ TT-) + P n< v - V-” • 

Here, the prime denotes deviation from the zonal mean. 

<K • K > is given explicitly by linearized equation as 
on 

oil 3V 

<W = - [ sr p n (u - u) + sr p n (u>v)] 

3U 3V„ 

- t7sl P n‘ u - v) + FSe P n< v ’ v >} 


3v v 


^ [U.P.(u.y) - V 0 P n (u,u)] 


r on 
31) 


av 


- C?T P n (u ’“> + 8F P n (v>u):l 


311 K 3V_K- 3w K 
r o n o n j. o n n 

- L ^TT* T T ^rv j > 


3x 


3y 3p 


where the suffix o denotes zonal mean. 


(4.21) 


(4.22) 
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V. RESULTS 

A. The Control Experiment 

In this section, we shall analyze the Kelvin waves, the mixed 
Rossby-gravity wave and the equatorially trapped Rossby wave, which 
are detected in the control experiment. 

1. Mean zonal wind 

Figure 2 shows the sunamer and winter mean zonal wind from the 
model simulation. The midlatitude tropospheric jets generated by the 
model are well positioned in both latitude and height compared with 
observations (Newell et al. , 1972) . The intensity of the jets is 
reasonably close to the observations, especially the Northern summer, 
while the jets in the Southern Hemisphere (SH) are 5 ms" 1 stronger 
than observed. A defect of the model is the weakness of the simulated 
tropical easterlies in both seasons. Instead of occupying the whole 
tropical troposphere, the tropical easterlies in the summer season are 
broken in the lower part of the troposphere and shift about 10° 
northward near the 200 mb level, while the surface easterlies and the 
easterlies near the tropopause compare favorably with observations. 

On the other hand, the tropical easterlies in the winter season 
disappear completely except for the weak surface easterlies simulated 
in agreement with observations. 

According to Maruyama and Yanai (1967) and Maruyama (1967) , mixed 
Rossby-gravity waves are generally observed when the prevailing flow 
is westerly, while Kelvin waves are observed in prevailing easterly 
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flow (Wallace and Kousky, 1968; Kousky and Wallace, 1971) . Thus, the 
weakness of the model tropical easterlies may have sane influences on 
the simulated equatorial waves. 

2. Kelvin wave in the summer season 

a. Spectral analysis Figure 3a shews the distribution of the 
power spectral density of the geopotential height at 2N and 200 mb 
during summer 1977 in the wavenumber-frequency domain. The power 
spectra demonstrate that EMWs (Eastward Moving Waves) are more 
pronounced than WMWs (Westward Moving Waves) . It is evident that 
transient wave energy is contained mainly in eastward propagating long 
waves. The most distinct spectral peak occurs at zonal wavenumber m=l 
and an eastward moving period in the range of 18-25 days. Another 
weaker spectral peak appears at m=2 and eastward period 5-6 days and 
has an energy content about a quarter of that of the strongest peak. 

Figures 3b and 3c show the power spectra of geopotential height 
and zonal wind for m=l as a function of frequency and latitude at the 
200 mb level. Distinct tropical and extratropical signals are 
distinguishable. The extratropical height variance in Figure 3b 
rapidly dominates beyond 10° latitude. Near the equator, the 
spectrum is dominated by an eastward propagating disturbance, which is 
symmetric with respect to the equator and decays away from the 
equator. Similar features can be also found in the spectrum of the 
zonal wind in Figure 2c. However, the center of the wind spectral 
peak is shifted southward and located at about 3S. Spectral peaks of 
both z and u occur at eastward periods of 18-25 days, which is 
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Figure 3. Power spectral analyses of (a) geopotential height (2N, 
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consistent with the spectral peak appearing at n*=l in Figure 2a. The 
spectra of meridional wind (not shown) are dominated by westward 
propagating disturbances in the extratropical region, while the 
spectral peak of the meridional wind corresponding to those of z and u 
is negligible. 

Figures 3d and 3e show height-frequency sections of the 
gecpotential height and zoned wind at 2N for m=l. Hie vertically 
coherent structure of eastward propagating disturbances with periods 
corresponding to that of the spectral peak identified above is 
evident. Careful observation of this disturbance indicates that the 
period of the disturbance shifts slightly towards shorter periods with 
increasing height. Hie maximum spectral peak of the geopotential 
height occurs at an eastward moving period of 20 days at 400 mb and 
approaches 15 days at 100 mb. In contrast to the power spectrum for 
z, the power spectrum for u for the EM* with a period of 18-25 days 
attains maximum at the 200 mb level. 

Features similar to those mentioned above can also be observed 
from the EMW associated with the spectral peak at ro=2 and periods of 
5-6 days (not shown). 

All the features that have been identified above are good 
indications of Kelvin waves with periods of 18-25 days for m=l, 
and with periods of 5-6 days for m=2. 

Since the two wave signals isolated from the spectral analysis 
have similiar behavior, we shall focus only on the signal at m=l 
and periods of 18-25 days in the following discussion. 
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b. Propagation characteristics Having isolated the wave 
signals, we are able to construct their horizontal structure as 
simulated by the model. Figure 4 presents the horizontal 
distributions of geopotential height and winds at 200 mb for m=l and 
periods of 18-25 days. Evidently, the simulated disturbance has a 
horizontal structure resembling that of the atmospheric Kelvin waves 
depicted by Matsuno (1966) , which are equatorially trapped and 
symmetric with respect to the equator in both pressure and zonal wind, 
with a negligible meridional wind component. Therefore, this 
disturbance is identified as a Kelvin wave with nt=l and periods of 18- 
25 days. 

Given the zonal wavenumber and Doppler-shifted frequency w , the 

latitudinal width Y L and vertical wavelength 1^ of the Kelvin wave may 
be computed as (Holton, 1979) 


2fi 

13 m 



(5.1) 


where S = 2. 2 8x1 0 -11 nf 1 s _1 , and S=4xl0 -4 s -2 is the square of the 

Brunt-Vaisala frequency. Estimated values of Y. and L for the Kelvin 

ij z 

waves are listed in Table 3. For example, the Kelvin wave with iw=l 
and period 18-25 days has a vertical wavelength of 5.9 km and a e- 
folding latitudinal width of 11° latitude, which is in agreenent with 
the wave pattern shown in Figure 3. 
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KELVIN WAVE (m=l) JUNE 28 



Figure 4. Horizontal distributions of geopotential height and 

horizontal winds at 200 mb for Kelvin wave with m=l and 
periods of 18-25 days in the control experiment oyer the 
northern summer. Contour interval is 1.5 m. Maximum 

wind at the equator is about 5 ms“* 
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Zangvil and Yanai (1980) and Yanai and Lu (1983) presented 
spectral studies of 200 mb level winds in the latitude belt 25S-45N 

for June- August 1967 and June-August 1972. They found Kelvin waves at 
ra=l with periods of 30 days and 7 days, and m=2 with periods of 7 days 
in the summer 1967. However, the wave signals in the summer 1972 are 
much weaker than those in the summer 1967 and only Kelvin waves at m=l 
with periods longer than 20 days and m=2 with periods near 15 days 
were observed. The great difference in equatorial wave activity in 
the two years is considered to be due to the interannual variation of 
the tropical general circulation during these two years. The year 
1967 had near normal monsoonal rainfall, while 1972 was a well-known 
drought and El Nino year. It was pointed out by Arkin (1982) that the 
tropical easterly jet is weaker and farther south in the summer of the 
year when the anomalously warm sea surface water associated with the 
Southern Oscillation appears over the equatorial Pacific (Warm Event) . 
Choi and van Loon (1987) further pointed out that the tropical 
divergent circulations, i.e. , the east-west walker and local Hadley 
circulations, during such summers are weakened and shifted eastward 
and that the changes in the tropical divergent circulations may cause 
a change in the energetic maintenance of the tropical easterly jet. 
Therefore, the change of wave activity between the two years 
suggests that the selection of the equatorial waves is related to 
global-scale forcing in the tropics. 

Yanai and Lu (1983) also calculated the equivalent depth and 
latitudinal width of the Kelvin wave which they had identified from 
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the observations. Comparing with their results, we find the Kelvin 
wave simulated by the CCM1 agree well with the observed ones as far as 
the horizontal and vertical structures of the waves are concerned. 

The simulated Kelvin waves are all confined within 16 ° latitude and 
have vertical wavelength less than 12 km, which suggests that these 
waves are generated in the upper troposphere instead of being free 
modes such as those proposed by Matsuno (1966) . In general, the 
observed Kelvin waves above the tropopause have larger vertical 
wavelengths than those of the upper tropospheric Kelvin waves (Wallace 
and Kousky, 1968; Hirota, 1978, 1979) . Salby and Garcia (1987) 
demonstrated that Kelvin waves with shorter periods in the equatorial 
atmosphere become relatively more important at upper levels since 
shorter periods are associated with larger vertical wavelengths and 
group velocities, and are less subject to radiative damping. 

Therefore, it is not surprising to have Kelvin waves with smaller 
vertical wavelengths and longer periods in the upper troposphere. 

Because of the distinct spectral character of the Kelvin wave 
signals it is possible to construct height-time cross sections for the 
time evolution of the Kelvin wave by band-filtering the eastward 
propagating features. Vertical structures of geopotential height and 
zonal velocity for m=l with periods of 18-25 days at 2N are shown in 
Figure 5. The oscillations of both z and u are clearly detectable and 
the vertical structures show a vertical phase change at the 500 mb 
level, which differs markedly from the structure of midlatitude 
systems characterized by strong baroclinicity. In the tropical 
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precipitation zone, convergence is not limited to low-level frictional 
inflow in the planetary boundary layer, but extends up to nearly 400 
mb level which is the height where the hot tcwers achieve their 
maximum buoyancy, while divergence attains its maximum at the 200 mb 
level. These features are clearly shown in Figures 5a and 5b. The 
disturbance has little phase tilt but a phase change somewhere between 
the 40Q mb and 500 mb over the entire sampling period. The maximum 
disturbance appears near 200 mb, where the maximum divergence occurs, 
implying the location of a possible energy source for Kelvin waves. 

There are two episodes of Kelvin wave disturbances during the 
sampling period. In the first episode, which covers days 30 to 80, 
the vertical structure of the disturbance is barotropic with a maximum 
amplitude of about 4 m at the 200 mb level. As time increases, a 
vertical phase shift appears near 400 mb while the amplitude increases 
to 8 m and the height disturbance extends up to the lower stratosphere 
with a slight phase tilt above 200 mb. On the other hand, the zonal 
wind component attains its maximum amplitude at 200 mb at the 
beginning of the episode and then the disturbance propagates downward. 

The second episode occurs during days 105-150. The disturbance 
of the geopotential height in this episode has a vertical phase change 
near the 500 mb level, which suggests that strong cumulus convection 
in the tropics occurs during this period. The amplitudes of z and u 
are about 8 m and 4 ms -1 , respectively. The fluctuations can be 
traced all the way up to the lower stratosphere. In the discussion of 
Figure 5, we also noted that the oscillations of the geopotential 
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height and the zonal wind are in phase. 

c. Precipitation and forcing inference Although the structure 
of Kelvin waves is now quite well understood, the nature of the 
tropospheric source for those waves is still somewhat uncertain. A 
number of possible tropospheric sources might be responsible for 
generating the observed waves. For example, geostrophic adjustment 
processes, nonlinear interaction among tropical tropospheric 
disturbances and forcing from middle latitudes all might play a role. 
However, the most likely energy source appears to be in-situ cumulus 
convection organized by synoptic-scale disturbances in the equatorial 
region. 

According to the satellite observations of cloud coverage 
compiled by Sadler (1970) and the satellite infrared effective 
temperature analyzed by Liebman and Hartmann (1981) , intense 
convection is associated with the intertropical convergence zone 
(ITCZ) , South Pacific Convergence Zone (SPCZ) , and the three tropical 
continents. Hie three regions are Central America and the northern 
part of South America, equatorial west Africa, and the monsoonal areas 
during the northern simmer. Chen (1985) pointed out that the 
stationary divergent modes, mainly the local Hadley and walker 
circulations, are responsible for the local maintenance of the high 
water vapor content over these tropical areas. Thus a high 
precipitation rate must exist over these areas due to the intense 
convection supported by the high water vapor content. 

Figures 6a and 6b present the model's mean-daily distribution of 
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the precipitation rate. Contour intervals are 2.5 nm day" 1 
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Figure 7. Latitude-time cross sections of (a) geopotential height 
for m=l at 200 mb and (b) precipitation rate for m=l-15 
with eastward moving periods of 18-25 days for the 
control experiment in simmer. Contour intervals are 2 m 
in (a) , and 3 nm day - * in (b) . Values greater than zero 


are hatched 
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precipitation and the root mean square of precipitation in the summer 
1977. Although the precipitation rate in the model (Figure 6a) 
generally tends to be somewhat larger than observed (i.e. , Pitcher et 
al. 1983) , the model successfully produces the general features of 
precipitation found in observations, except for the precipitation 
deficit over the equatorial Africa. For example, the intense 
equatorial rain belts, including an elongated maximum over India and 
the western Pacific in accordance with the Asian monsoon circulation, 
and a tongue of rainfall in the Central America, extending from the 
northeast Pacific to the north part of Central America, are reproduced 

in the model. It is evident that the daily rainfall fluctuates 

\ 

substantially over the three maximum precipitation centers (Figure 
6b). 

Figure 7a shows the latitude-time cross section of the 200 mb 
geopotential height for the Kelvin wave and the precipitation 
disturbances with periods of the Kelvin wave. Hie time evalution of 
the Kelvin wave associated with m=l and periods of 18-25 day at 150E 
longitude can be seen clearly in Figure 7a. It is evident that the 
wave is confined within an equatorial belt between 10S to ION during 
the entire sample period. Hie two episodes of the Kelvin wave, which 
were identified before, can also be observed. 

As pointed out by Holton (1972) , there is no convincing evidence 
to show that the field of di aba tic heating in the tropical troposphere 
has a spectral peak at the frequency where Kelvin waves are observed. 
Hie spectral analysis of the model's precipitation rate (not shown) 
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gives a stochastic distribution in the wavenumber-f requency domain. 
Based on this argument, Figure 7b shows the precipitation perturbation 
associated with m=l-15 and eastward periods of 18-25 days. The 
interesting feature of the rainfall disturbance is the northward 
migration of its maximum center. For example, at day 26, the maximum 
center is located near 10S. As time progresses, it continuously moves 
northward until it approaches the equator at day 105. Recalling the 
vertical structures of the Kelvin wave shown in Figure 5, we find that 
the vertical phase change of the geopotential height after day 105 
coincides with the precipitation disturbance near the equator. This 
finding strongly suggests that the atmosphere is more sensitive to 
forcing in the frequency band of Kelvin waves. 

In order to understand the possible sources of excitation for 
wave disturbances in the tropics, Holton (1972, 1973) and 
Murakami (1972) studied the atmospheric response to prescribed 
tropospheric heat sources in linearized primitive equation models. 
Holton (1972) simulated the atmospheric response to the latent heat 
released in the ITCZ. In any linearized model, the atmospheric 
response to forcing in one hemisphere is equivalent to the sum of the 
responses that are symmetric and antisymmetric with respect to the 
equator. Therefore, Holton studied the atmospheric responses of 
symmetric and antisymmetric heating separately. He concluded that a 
heating oscillation which is symmetric about the equator will generate 
eastward moving Kelvin wave and the frequencies of the Kelvin waves 
are determined by the natural bandpass selectivity of the atmosphere 
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combined with a red noise forcing spectrum in time or space. 

On the other hand, Murakami (1972) simulated a quite realistic 
Kelvin wave as the response to an eastward moving heat source with 
m=l-2 and a period of about 16 days at the latitude of the monsoon 
region (20N) . This specification of heat source is partly supported 
by observational evidence, which shows that during the monsoon season 
the rainfall fluctuates substantially at that latitude with periods 
longer than 10 days, and a major portion of zonal rainfall asymmetry 
is presumably carried by wavenumbers 1 and 2. 

As revealed in Figure 7, it does not matter whether the eastward 
propagating heat source is located south of the equator or straddles 
the equator; the Kelvin wave is observed in both cases. This finding 
is in agreement with both Hoi ten's and Murakami's results. 

An attempt was next made to determine the effect of a zonally 
asymmetric heat source on the generation of Kelvin waves. Figure 8 
shows the amplitude of the Kelvin wave in geopotential height at 200 
mb and the preciptation rate with im*l-15 and the period of the Kelvin 
wave at (90W, 2N) , (105E, 2N) and (150E, 2N) , respectively. The solid 
lines indicate geopotential height and the dashed lines precipitation 
rate. Hie Kelvin wave disturbance and precipitation in the eastern 
equatorial Pacific (90W, 2N) or in eastern Indian Ocean (105E, 2N) 
show little correlation. However, the two fields show a high 
correlation in the western Pacific Ocean (105E, 2N) . 

We also noted that there is a time lag between the geopotential 
height and precipitation oscillations at a given location. This time 
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lag singly implies an eastward tilt of the wave structure with 
increasing height. From theoretical considerations, the axes of 
maximum upward motion and maximum zonal wind must be in phase in order 
to satisfy the continuity equation. Consequently, precipitation 
caused by vertical motion is also in phase with zonal wind 
oscillations. On the other hand, since the zonal wind component is 
geostrophically balanced by the meridional pressure gradient, it is 
expected that the maximum zoned winds and maximum pressure should 
coincide. Thus, at a given longitude, the pressure oscillation at 
upper levels will lead the lower level precipitation oscillation as 
long as the vertical structure of the Kelvin wave tilts eastward. 
Therefore, Figure 8 strongly suggests that the maritime region of 
Indonesia in the western Pacific may be a possible region where Kelvin 
waves originate. 

d. Maintenance of the Kelvin wave Since the data output from 
CCM1 is interpolated to conventional meteorological levels, the 
maintenance of a particular wave disturbance can be gained by 
computing every term in the kinetic energy equation (4.17) , which is 
proposed by Hayashi (1980) , in the wavenumber-frequency domain. In 
general, the meridional wave energy flux due to the pressure work and 
the poleward momentum transport, which is included in the nonlinear 
term, have maxima at 200 mb, simply because this level has the maximum 
kinetic energy content of wave motion. 

The first term on the right-hand side of the equation (4.17) 
represents transfer of kinetic energy from or to wavenumber n through 
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nonlinear dynamic interaction. Hie second term represents the 
conversion of eddy available potential energy to eddy kinetic energy. 
Hie third and fourth terms represent the meridional and vertical 
convergence of the potential energy flux associated with wave 
disturbances. 

Figure 9 presents the energetic budget analysis of equation 
(4.17) as a function of latitude and frequency for m=l. It reveals 
that the kinetic energy (Figure 9a) associated with the Kelvin wave 
with a centered period of 20 days is dominant. However, the 
latitudinal structure of the kinetic energy is not symmetric with 
respect to the equator. Hie maximum kinetic energy is located at 2S. 
It suggests that the asyirmetry of the basic zonal flow (Figure 2) may 
be responsible for the asymmetric distribution of Kelvin-wave kinetic 
energy. 

Hie conversion of eddy available potential energy to eddy kinetic 
energy is shown in Figure 9b. TVo significant energy conversion peaks 
appear near 12-day and 6-day periods. It is seen that a positive 
energy conversion with a period about 12 days is located south of 5S 
and a negative energy conversion with the same period appears north of 
5S. Using station data of upper air observation over the Marshall 
Islands area (165E, 5N) from April 15 to July 31, 1956, Nitta (1970) 
computed the vertical distribution of conversion of available 
potential energy into kinetic energy and found that a large energy 
conversion from eddy available potential energy to eddy kinetic energy 
is concentrated in two period ranges centered at 12 days and 5 days 
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Figure 9. Various energy conversions (itfI, eastward moving 

periods) at 200 mb for the control experiment over the 

9P, (<)> >w) 

northern summer: (a) K. , (b)- P, (a,w ) , (c) L , 

3P (a v ) 1 ?P 

(d)- and (e) N.(k) . Contour intervals are 2 
9y 1 

J in (a); 0.5xl0~ 5 W in (b) , (c) and (e); and 10" 5 W in 
(d) . Values greater than zero are hatched 
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with the maximum energy conversion being located at the levels between 
400 mb and 200 mb. Thus, he suggested that the conversion of 
available potential energy is the most important one for the 
maintenance of large-scale disturbances in the tropics. However, 
Nitta's argument does not seem to be the case for the Kelvin wave 
since the period of latter is longer than that of Nitta's waves. 

Figure 9c shows the vertical convergence of potential energy 
flux. It is evident that the eastward component of this quantity is 
positive (convergence) for almost the entire period at this level. 

The maximum convergence occurs at 2N with a period range centered at 
20 days, which coincides with the period of the Kelvin wave. Hayashi 
et al. (1984) showed that the eastward moving component of the 
vertical potential energy flux transports upward above 7km while below 
that level it transports downward. The reversal of the direction of 
energy flux at 7 km suggests that the wave energy source exists near 
this level. In Figure 9c, the most distinct convergence of potential 
energy flux appears north of the equator. From the discussion in the 
previous section, we can infer that the local divergent circulations, 
which are generated by deep penetrative convections in the tropics, 
are responsible for this vertical energy transport. Thus, cumulus 
convection in the tropics is important for the maintenance of the 
Kelvin wave. 

Figure 9d shows the meridional convergence of energy flux. The 
negative area between 8S to 5N with periods 15-30 days indicates 
divergence of the meridional energy transport, while the positive area 
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poleward of the negative area implies convergence of the meridional 
energy flux. This result confirms the finding obtained by Yanai and 
Lu (1983) , i.e. , that there is a meridional divergence of energy flux 
associated with the Kelvin wave. Fran Figures 9c and 9d we may 
conclude that the energy source at the 200 mb level is primarily the 
convergence of vertical energy transport to the north of the equator, 
ftie net energy input is then transported away from the equator to 
maintain the Kelvin wave. 

In order to understand the basic dynamics of equatorial waves , 
most of the previous studies of this subject were based upon the 
linear system. Since real atmospheric motions are not linear, it is 
of interest to see the importance of nonlinear interaction in the wave 
maintenance. Figure 9e shows the transfer of kinetic evergy into EMtfs 
with m=l through nonlinear interaction. It is seen that the 0W 
associated with m=l and periods of 12-30 days gains kinetic energy 
through nonlinear processer in the SH, while the wave loses kinetic 
energy in the Northern Hemisphere (NH) between the equator and ION. 

The nonlinear energy transfer can be further partitioned into 
wave-mean flow and wave-wave interactions, which are shown in Figures 
10a and 10b, respectively. The transfer of kinetic energy through the 
wave-mean flow interaction is significant for those EMVs with periods 
longer than 12 days. It is evident that EMfe with periods longer than 
18 days, in which the Kelvin wave is observed, gain kinetic energy 
from zonal flow within a latitudinal belt between 5S to 5N (Figure 
10a) and loses kinetic energy in the NH through the wave-wave 
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interaction (Figure 10b) . In addition, the EJW with a period of 15 
days gains kinetic energy in the SH through the wave-wave interaction 
and then transfer the energy into zonal flow through the wave-mean 
flow interaction. 

3_. Kelvin wave in the winter season 

In this section we shall briefly discuss the characteristics of 
the Kelvin wave simulated by the model in the winter season with 
emphasis on the differences between the Kelvin waves in the summer 
and those in the winter. 

a. Spectral analysis The power spectrum of z at 2S and 200 
mb is shown in Figure 11a as a function of wavenumber and frequency. 

It is evident that the power spectrum of z for EMWs has a spectral 
peak at m=l and periods of about 15-20 days. A weak spectral peak 
associated with WMWs occurs at m=3 and periods of 20-30 days. This 
signal is later identified as an equatorial trapped Rossby wave. 

Figures lib and 11c show the meridional structures of both z and 
u with m=l. In Figure lib, both WMtf and EWV are dominant in the 
extratropical area. A weak spectral peak is confined in the tropics 
between 8S-5N with eastward periods of 15-20 days. This spectral peak 
corresponds to the one at m=l and periods of 15-20 days in Figure 11a. 
The power spectrum of u with itF=l (Figure lib) has larger power 
concentrated in the low frequency part of the spectrum, while 
appreciable power of EMWs appears near the equator with periods of 15- 
20 days. 
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Figures 5d and 5e show the vertical structures of P(z) and P(u) 
associated with n*=l at 2S. Both eastward variances of z and u are 
significant at the period range of 15-20 days. Appreciable powers are 
evident from 300 mb up to 100 mb and from 850 mb up to 100 mb in 
Figures lid and lie, respectively. 

b. Propagation characteristics The power spectrum of 
wavenumber 1 exhibits distinct spectral features corresponding to 
eastward propagation with periods of 15-20 days, consistently over a 
band of tropical latitudes and several pressure levels. Figure 12 
shows that the disturbance at 200 mb has a horizontal structure 
strongly resembling that of atmospheric Kelvin waves. From linear 
wave theory, the vertical wavelength and latitudinal width of the wave 
are estimated to be 6.3 km and 12° latitude, respectively (Table 3). 
Comparing with the Kelvin wave in the summer season, we find the 
propagation characteristics of the Kelvin waves in both seasons are 
essentially the same, except for the difference in periods associated 
with the waves, which might be due to the seasonal change of the 
heating field organized by large-scale motions. 

Tsay (1974) examined the characteristics of the equatorial waves 
identified in a perpetual January run with NCAR global circulation 
model. Hie Kelvin waves simulated by the model were found at m=l with 
an eastward traveling period of about 24 days and m=2 with an eastward 
period of about 12 days. Hie vertical wavelength and latitudinal 
width of the Kelvin wave with m=l were 8 km and 14° latitude, which 
agreed reasonably well with those of observed waves (Wallace and 
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Kousky, 1968? Kousky and Wallace, 1971) . However, the observed 
oscillation with m=l has a significantly shorter period (15 days) . 
Furthermore, no response at wavenunber 2 at one-half* the wavenumber 1 
period was detected in the observations. On the other hand, the 
Kelvin wave simulated by the 0CM1 is even closer to the observed wave, 
comparing to Tsay's study, as far as the latitudinal width, period and 
vertical wavelength are concerned. 

The forcing of the Kelvin wave in the winter season can be 
inferred from Figures A1-A4. 

Cm Maintenance of the Kelvin wave Figure 13 shows the 
maintenance of the Kelvin wave at 200 mb in the winter season. The 
most distinct feature is that the conversion of available potential 
energy into kinetic evergy plays an inportant role in maintaining 
Kelvin wave (Figure 13b) . The maximum conversion occurs at periods of 
about 15-20 days and near 3 S, while the convergence of the vertical 
energy flux exhibits less importance in the wave maintenance (Figure 
13c) . This result agrees with the conclusion obtained by Manabe et 
al. (1970) that the conversion of available potential energy generated 
by the heat of condensation is most inportant for the maintenance of 
large-scale disturbances in the tropics. 

The meridional wave energy flux associated with the Kelvin wave 
is characterized by convergence in the SH and divergence in the NH 
(Figure 13d) . The former probably results from the conversion of 
available potential energy into kinetic energy near 3S. This outcome 
is contrary to the conclusion obtained by Tsay (1974) who found that 
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the energy source of Kelvin wave in the layer 12-18 km is due 
primarily to convergence of equatorward energy transport. 

Nonlinear kinetic energy transfer associated with EMWs at m=l is 
show in Figure lie. It is evident that a large energy transfer take 
place in the period range of 12-30 days, in which the Kelvin wave is 
observed. In contrast to the simmer case, there is positive energy 
transfer in the NH and negative energy transfer in the SH. 

We have seen that the maintenance of the Kelvin wave in the 
winter is somewhat different from that in the summer. The difference 
might be caused by the seasonal change in the heating field and the 
mean zonal wind distribution. The former may alter the global 
divergent circulation and the latter affects the wave-mean flew 
interaction. 

4. Mixed Rossby-gravity wave in the summer season 

a. Spectral analysis The spectrum of the meridional wind 
component at 2N and 200 mb as a function of wavenumber and frequency 
is shown in Figure 14a. There are severed spectral peaks at m=5 , among 
which the most distinct spectral peak is associated with westward 
moving periods of 12-15 days and characterized by a unique spectral 
peak in the zonal wind component at the equator (Figure 3a) . Shown in 
Figures 14b and 14c are the latitudinal distributions of P(z) and P(v) 
associated with m=5. P(z) is dominated by both EMWs and WMWs in 
extratropical area, while the WMWs at periods of 12-30 days show a 
very clear equatorward decay. On the other hand, P(v) attains a 
maximum near the equator, especially for the WMWs with periods of 12- 
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Figure 14. Power spectral analyses of (a) meridional wind (2N, 200 
irb) , (b) wavenumber 5 geopotential height (200 mb) , (c) 
wavenumber 5 meridional wind, (d) wavenumber 5 
geopotential height (15N) and (e) wavenumber 5 meridional 
wind (2N) for the control experiment over the northern 
summer. Contour intervals are 1 irrs“2 day j. n ( a ) / 15 
it r day in (b), 1.5 nrs - ^ day in (c) , 8 nr day in (d) 
and 0.8 irrs “2 day in (e) 
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15 days (Figure 14c) . Figures 14d and 14e demonstrate the vertical 
distributions of P(z) and P(v) associated with m=5 at 15N and 2N, 
respectively. It is evident that almost all the spectral peaks occur 
near 200 mb. 

Knowing the horizontal and vertical structures of these spectral 
peaks, we infer that the the spectral peaks with eastward periods are 
eastward moving inertial-gravity waves which might be generated by the 
model. The nature of the spectral peak associated with westward 
periods of 5.5 days is not clear since the horizontal structure of the 
disturbance does not fit into any wave structure predicted by the 
linear wave theory. All the features associated with the spectral 
peak at m=5 and periods of 12-15 days are indications of a mixed 
Rossby-gravity wave. 

b. Propagation characteristics The horizontal distributions 
of geopotential height and horizontal wind components associated with 
m=5 and periods of 12-15 days at 200 mb are shown in Figure 15. The 
geopotential height distribution exhibits approximate antisymmetry 
about the equator with high and low centers located around 13N and 6S. 
Hie horizontal wind is clockwise around the pressure system of the NH 
low and the SH high. Hie estimated amplitude of the meridional wind 
at the equator is about 0.6 ms~' L . Hie amplitude of the mixed Rossby- 
gravity wave observed by Zangvil and Yanai (1980) is about 0.8 ms 
Yanai and Murakami (1970) obtained an amplitude for meridional 
velocity of about 2.5 ms“^ (this includes contributions from all 
wavenumbers of EMW and WJftf) at 200 mb. Thus the amplitude of the 
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Figure 15. Same as Figure 4, except for the mixed Rossby-gravity 
wave with m=5 and periods of 12-15 days 
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roixed Rossby-gravity wave in the present study is very close to those 
observed. 

For mixed Rossby-gravity waves, the latitudinal width Y L and the 
vertical wavelength L z may be computed fran (Holton, 1979) 
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Hie estimated values of Y^ and 1^ are given in Table 3. For example, 
from Figure 2, if we take the Doppler-shifted phase speed to be 5 ms" 1 
in the SH, then the latitudinal width is about 8° latitude and the 
vertical wavelength is about 3 km. 

Figure 16 shews a height-time section of meridional wind with m=5 
and periods of 12-15 days at 2N during the period in which the mixed 
Rossby-gravity wave is significant. It reveals a regular ocsillation 
in v associated with the mixed Rossby-gravity wave and exhibits a 
little tendency of downward propagating phase in the upper 
troposphere, which is in agreement with the result obtained by Tsay 
(1974) who showed that the mixed Rossby-gravity wave (ro=2, period=4.6 
days) simulated by the NCAR stratospheric GCM have a downward 
propagating phase. 
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Figure 16. Height-time cross section of meridional wind with m=l 
and westward moving periods of 12-15 days for the 
control experiment over the northern summer. Contour 
interval is 0.5 ms - *. Values greater than zero are 


hatched 
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c « Precipitation and forcing inference As mentioned in 
Chapter II, there are two iirportant hypotheses on the generation of 
mixed Rossby-gravity waves. One is the lateral forcing theory 
proposed by Mak (1969) , who assumed that the vertical wavelength of 
the wave coincides with the vertical wavelength of lateral forcing 
through linear resonance. 'Urns, the westward moving resonant waves 
have a large response at m=3 and periods of 4 days. However, Hayashi 
(1976) noted that by imposing midlatitude forcing as a lateral 
boundary forcing, the resonant waves are spurious "westward Kelvin 
waves" having a horizontal structure completely different from mixed 
Rossby-gravity waves. 

The other hypothesis is called wave-CISK theory. This theory 
proposes that the interaction between equatorial waves and cixnulus 
convection is unstable and that it can produce mixed Rossby-gravity 
waves. However, this theory predicts sharp spectral peaks in the 
heating field, which is not detected in the real atmosphere. Very 
recently, however, it has been shown by two studies (Hayashi and Sumi, 
1986; Lau and Peng, 1987) that if a nonlinear effect is taken into 
account in wave-CISK, results concerning the wavenumber and period 
selection change completely. Therefore, it would be interesting to 
investigate whether a nonlinear model with wave-cumulus interaction 
produces mixed Rossby-gravity waves with the same spectral peak in 
wavenumber and frequency as observed waves. 

Figure 17 shows latitude-time sections of the oscillations of z 
associated with the mixed Rossby-gravity wave and precipitation with 
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(a) Z (150E) 



(b) PRECIPITATION (150E) 
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Figure 17. Latitude-time cross sections of gecpotential height 

with m=l at 200 mb (a) and precipitation rate with m=l- 
15 (b) associated with westward moving periods of 12-15 
days for the control experiment over the northern 
summer. Contour intervals are 1 meter in (a) , and 2 mm 
day”-*- in (b) . Values greater than zero are hatched 
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m=l-15 and periods of 12-15 days. Figure 17a reveals clearly the 
oscillation of the mixed Rossby-gravity wave with the maximum pressure 
at 6S and 12N, especially during the period of days 45-57 (May 14-26) . 
The precipitation (Figure 17b) exhibits an antisymmetric structure 
with respect to the equator, while the maximum precipitation in both 
hemispheres migrates equatorward. From the two figures, it is evident 
that the oscillations of the two fields are in phase during days 45- 
57. 

The significant coherence between the mixed Rossby-gravity wave 
and precipitation indicates that latent heat release from cumulus 
convection might be the energy source for the mixed Rossby-gravity 
wave. This confirms the results obtained by Itch and Ghil (1988) , who 
recently performed a numerical experiment to investigate the 
excitation mechanism of the mixed Rossby-gravity waves in the tropical 
troposphere. They found the coherence between lateral forcing and 
low-level convergence, on one hand, and between the mixed Rossby- 
gravity wave and convergence, on the other, are very high. Therefore, 
the lateral forcing, highly correlated with low-level convergence and 
precipitation, appears to influence mixed Rossby-gravity waves via 
low-level convergence and ensuing precipitation, but exhibits only a 
low direct correlation with the waves. 

d. Maintenance of the mixed Rossby-gravity wave Figure 18 
presents the energetic budget analysis of equation (4.17) for the 
mixed Rossby-gravity wave. It is evident that the WMWs with m=5 have 
appreciable kinetic energy in the period range of 5-60 days. The 
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Figure 18. Same as Figure 9, except for wavenumber 5 and westward 
moving periods. Contour intervals are 0.8 J in (a); 
0.5xl0 -5 W in (b), (c) and (e); and 10“ 5 W in (d) . 
Values greater than zero are hatched 
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Figure 19. 
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Same as Figure 10, except for wavenumber 5 and westward 
moving periods 
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kinetic energy associated with the mixed Rossby-gravity wave is almost 
constant from 15S to 15N. 

The conversion of available potential energy into kinetic energy 
is displayed in Figure 18b. Although this quantity becomes 
significant toward the low frequency regime, the contribution from 
this quantity to the mixed Rossby-gravity wave over the period of 12- 
15 days can be almost neglected. This indicates that the wave has a 
somewhat barotropic nature as shown in Figure 16. 

Figure 18c shows the vertical convergence of wave energy flux. 
There are two convergence centers with periods of 12-30 days lying at 
12S and 6N, which imply two ascending branches of the Hadley cells 
associated with cumulus convection during Northern Summer. Again, it 
demonstrates that the mixed Rossby-gravity wave is maintained, to seme 
extent, by latent heat release in the tropics (Hayashi and Golder, 

1978) . 

The meridional convergence of wave energy flux shown in Figure 
18d is characterized by two convergence centers near the equator at 
periods of 5-6 days and 12-30 days, respectively. In agreement with 
previous studies by Zangvil and Yanai (1980) and Yanai and Lu (1983) , 
there is appreciable equatorward wave energy flux associated with the 
WMW at period range of 12-30 days, in which the mixed Rossby-gravity 
wave is observed. Although the resonant wave theory of Mak (1969) 
proved to be unrealistic, the previous studies mentioned above still 
show that the mixed Rossby-gravity waves are closely related to the 
energy flux from midlatitudes. 
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Figure 18e presents the kinetic energy transfer into the WMWs 
with ro=5 through nonlinear processes. It is evident that the 
nonlinear effects become significant for these waves with periods 
longer than 5 days. The nonlinear energy transfer associated with the 
mixed Rossby-gravity wave is positive north of 3N and negative south 
of 3N. 

Figure 19 shows the nonlinear transfer due to the wave-mean and 
wave-wave interactions. In Figure 19a, it is evident that the WMWs at 
periods of 12-15 days gain kinetic energy from the mean flow north of 
2N and south of 4S. However, based upon linear wave theory, Zangvil 
and Yanai (1980) showed that the mixed Rossby-gravity wave with m=4 
and a period of 5 days is accompanied by a conversion of eddy kinetic 
energy to kinetic energy of the mean flow. 

Ohe energy transfer due to the wave-wave interaction is displayed 
in Figure 1%). The energy transfer for those WMtfs with periods longer 
than 5 days is significant, while the WMWs with periods of 12-15 days, 
for which the mixed Rossby-gravity wave is observed, lose kinetic 
energy between 15S and 5N. 

5. Rossby wave in the winter season 

a. Spectral analysis Figure 20 presents the power spectral 
analysis of the Rossby wave. The top figure shows the power spectrum 
of z at the 200 mb level and 15N in the wavenumber-frequency domain. 
There is a spectral peak at m=3 and period range of 20-30 days. A 
similar pattern is also exhibited in P(u) at 2N (not shown). The 
vertical structure of P(z) with m=3 at 15N reveals a significant peak 
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Figure 20. Power spectral analyses of (a) geopotential height 
(15N, 200 mb ) , (b) wavenumber 3 geopotential height 
(15N) , and (c) wavenumber 3 zonal wind (2N) for the 
contrpl experiment over northern winter. Contour 

intervals are 25 m 2 day in (a) and (b) , and 1.5 n^s -2 
day in (c) 
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Figure 21. Same as Figure A, except for the Rossby wave with m=3 
and westward moving periods of 20-30 days in the 
control experiment over the northern winter. Contour 
interval is 2 m 
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with westward periods of 20-30 days at 200 mb. P(z) at this period is 
significant. A weaker spectral peak with eastward periods of 20-30 
days is also observed. The vertical distribution of P(u) with m=3 at 
2N behaves the same way as P(z) . A significant spectral peak with 
westward periods of 20-30 days occurs from 400 mb all the way up to 100 
mb/ while another spectral peak corresponding to the one with eastward 
period in P(z) is also noticed. The spectral analysis indicates that 
P(z) attains its maximum somewhere away from the equator, and that P(u) 
has its maximum on the equator. It strongly suggests that the spectral 
peak with westward period is a signal of the existence of a Rossby wave 
with m=3 and periods of 20-30 days. The nature of the spectral peak 
associated with eastward period is not clear. 

b. Propagation characteristics Figure 21 shows the 

■ — I ■■ ! — ■■ ■ 

latitudinal structures of geopotential height and horizontal wind of 
the Rossby wave for m=l and periods of 20-30 days. It is evident that 
z and u have symmetric structures and v has antisymmetric structure 
with respect to the equator, which corresponds to the Rossby wave with 
latitudinal nodal wavenumber n=l depicted by Matsuno (1966) . 

Following the correction method used by Diky and Golitsyn (1968) 
and Walkerscheid (1980) , we calculate the vertical wavelength and 
latitudinal width by taking into account the effects of mean zonal 
flow. From Figure 2, if we take the mean zonal flow to be 10 ms - ^ in 
the tropics, then the effective rotation rate of the earth is &e = 

1.02 Oo, D o is the rotation rate of the earth and the Doppler-shifted 
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angular frequency of the Rossby wave w =0.11rt o . Thus the vertical 
wavelength of the Rossby wave estimated from the dispersion curves of 
Longuet-Higgins (1968) is about 16 km and the maxi mm pressure is 
located at 15° latitude. Comparing with the Figure 21/ we find that 
the location of the maximum pressure in the SH is consistent with that 
obtained from linear wave equation, while the one in the NH is not 
well predicted theoretically. This wave is identified as an 
equatorially trapped Rossby wave with zonal wavenumber m=3, 
latitudinal wavenumber n=l and a period of 20-30 days. 

Although Rossby-type waves in the upper troposphere have been 
noted by synoptic analysis (Yanai, 1963) and by spectrum analysis 
(e.g. , Yanai and Murakami, 1970; Nitta, 1972), the evidence from 
time-spectrum analysis was not strong. Gruber (1975) presented seme 
evidence suggesting the existence of a Rossby wave with m=5-6 and a 
period of about 7 days. Yanai and Lu (1983) identified an 
equatorially trapped Rossby wave with m=2 and a period of 10 days, and 
with m=2 and a period of 30 days from observational data at 200 mb 
during 1967 and 1972, respectively. Hie amplitudes in u of both waves 
are about 3-4 ms - ^ at the equator, which is close to 2.7 ms“^ of the 
present study. 

Figure 22 shews the height-time section of the oscillation in z 
associated with the Rossby wave at 150E longitude and 15N. It reveals 
a regular oscillation of the Rossby wave with maximum amplitude near 
the 200 mb level. 
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c. Maintenance of the Rossby wave Figure 23 presents the 
energetic budget analysis of equation (4.17) for the Rossby wave. The 
kinetic energy with m=3 is shown in Figure 23a. Large kinetic energy 
is concentrated in the long-period part of the spectrum. This kinetic 
energy with periods of 20-30 days, at which the Rossby wave is 
observed, is significant from 3 OS to 3 ON with a maximum at ION. 
Midlatitude system s with a period of about 30 days are also observed 
at 3 OS and 30N. 

Figure 23b shows the conversion of available potential energy 
into kinetic energy. Significant conversions take place in two period 
ranges: 5-7 days and 15-30 days. We note that the period range of 15- 

30 days corresponds to that of the Rossby wave. The vertical 
convergence of wave energy flux in this period range is also 
significant (Figure 23c) . Over this period range, there are two 
convergence centers lying at 20S and 10N, where ascending branches of 
winter Hadley cells are located. Therefore the thermal effect is one 
of the important factor for the maintenance of equatorially trapped 
Rossby waves. 

Figure 23d presents the meridional convergence of wave energy 
flux. It is evident that the convergence associated with WMWs with 
periods less than 5 days is very weak. However, the meridional 
convergence of wave energy flux for those WMWs with periods of 20-30 
days, including the Rossby wave, is significant between the equator 
and 20S, and north of 15N. Between the two regions of convergence, 
wave energy flux diverges at 10N transporting out the energy 
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Figure 23. Same as Figure 9 , except for wavenumber 3 and westward 
moving periods. Contour intervals are 1.5 J in (a); 
0.5xl0" 5 W in (b); 10" 5 W in (c) , (d) and (e) 
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Figure 24. Same as Figure 10, except for wavenumber 3 and westward 
moving periods. Contour intervals are 0.5x10“^ W in 
(a) , and 10“^ w in (b) 
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accumulated by the thermal field at this latitude. We also notice 
that a divergence of wave energy flux with periods of 12-30 days is 
located at 30S. This divergence indicates disturbances from the 
midlatitude of the SH can freely propagate into the tropics as 
discussed by Charney (1969) and Bennett and Young (1971) . 

The kinetic energy transfer into WMWs with m=3 through nonlinear 
processes is displayed in Figure 23e. Similarly, appreciable energy 
transfer takes place only in the period range of 12-30 days, in which 
the Rossby wave is identified. It is evident that the Rossby wave 
loses kinetic energy at all the latitudes except near ION, where the 
wave gains energy through nonlinear processes. 

Separating the nonlinear energy transfer into wave-mean flow and 
wave-wave interactions (Figure 24) , we find that kinetic energy is 
transferred into the Rossby wave near ION in both nonlinear 
interactions, while the Rossby wave gains kinetic energy between 20S 
and 5S in the wave-mean flow interaction and then the kinetic energy 
is transferred into other waves in the wave-wave interaction. 

B. The No-mountain Experiment 

From the above analysis of the control experiment, we have seen 
that the NCAR CCM1 is able to simulate equatorial planetary waves 
which compare reasonably well to observations. In order to evaluate 
the role of the earth's orography in the excitation of these waves, we 
analyzed a numerical experiment with CCM1 by eliminating mountains 
from the model. As the result, only Kelvin waves are identified in 
both winter and summer seasons. Since the Kelvin waves in both 
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seasons have similar propagation characteristics , we shall only 
discuss the Kelvin wave in the summer season in this section. 

1_. Spectral analysis 

Figure 25a shows the power spectrum of z at 200 mb and 2N in the 
wavenumber-frequency domain. A striking feature of this figure is the 
concentration of P(z) in the low-wavenumber and long-period regime. 
Examining the effect of the equatorial Pacific sea surface temperature 
anomaly on the atmospheric circulation using the NCAR CCMO with and 
without mountains, Blackmon et al. (1987) found that removing the 
mountains in the model produces a more zonally uniform flow and the 
stationary waves have smaller anplitudes than those in the mountain 
case, especially for zonal wavenumbers 2 and 3. Conparing with Figure 
25a and Figure 3a, we find the power of z associated with wavenumber 2 
and 3 is reduced significantly. It is evident that EMNs are dominant, 
especially the disturbance with m=l and periods of 15-20 days. 

Shown in Figures 25b and 25c are P(z) and P(u) with m=l at 200 mb 
as a function of latitude and frequency. Both spectra are dominated 
by EMWS, which are symmetric with respect to the equator and decay 
away from the equator, while the power spectrum of the meridional wind 
at m=l is extremely small (not shown) . We also note that the 
extratropical disturbances, especially stationary disturbances, are 
reduced markedly in the no-mountain experiment. Therefore, the 
possible interaction between midlatitude and tropics may be reduced. 

Figures 25d and 25e show height-frequency sections of zonal wave 
1 geopotential height and zonal wind component at 2N. The vertically 
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coherent structures of the EMWs at the period identified above are 
evident. 

2. Propagation characteristics 

The above spectral analysis strongly suggests the possible 
existence of the Kelvin wave at m=l and periods of 15-20 days. Figure 
26 shows the horizontal distributions of geopotential height and winds 
at 200 mb for periods of 15-20 days with m=l. From equations (5.1) 
and (5.2) , the latitudinal width and vertical wavelength in this case 
are estimated to be 13° latitude and 7 km, respectively. In Figure 
26, the amplitudes of fluctuations of zonal wind are 1-3 ms -1 , and the 
fluctuations of z are 10-14 m at the equator. It is evident that the 
characteristics of the Kelvin wave in this case are similar to those 
in the control case except with a smaller period. This suggests that 
the topography, with which the stationary wave and midlatitude 
disturbances are associated, has little effect on the propagation 
characteristics of Kelvin waves in the tropics. 

Figure 27 shows the oscillations of the Kelvin wave in z and u at 
2N. It reveals that there are two episodes of the Kelvin wave in this 
sample period: days 0-50 and 95-130. In contrast to those in the 

control run, the eastward propagating modes in this case exhibit 
downward phase progression over the entire sample period, with the 
maximum amplitudes occurring in the lower stratosphere. Hie two 
features mentioned above indicate the effect of topography on large- 
scale motions. As pointed out by Held (1983) , thermally forced eddies 
and mountain forced eddies have rather different structures. Hie 
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Figure 26. Same as Figure 4, except for the Kelvin wave with m=l 
and periods of 15-20 days for the no-mountain 
experiment. Contour interval is 2 m 
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former are of larger scale than the latter, particularly in the upper 
troposphere. They also have a larger and more systematic westward 

i 

phase shift with height and show somewhat less of a tendency toward 
maximum amplitude at the tropopause. The mountain forced eddies, on 
the other hand, tend to reach their maximum amplitudes near 200 mb and 
exhibit less westward phase tilt with height. 

Vertical phase change of the zonal wind oscillation during the 
episodes of the Kelvin wave (Figure 27b) indicates that the vertical 
structure of low-level convergence and upper-level divergence 
accompanies the Kelvin wave. Corresponding features can also be 
observed in Figure 27a. 

3. Precipitation and forcing inference 

Mountains have a considerable effect on the longitudinal 
distribution of tropical precipitation. Figures 28a and 28b show time 
averaged precipitation rate and the precipitation variance for the no- 
mountain case. Removing mountains from the model produces a more 
zonally uniform distribution of tropical rainfall, which agrees with 
the results obtained by Manabe and Terpstra (1974) and Blackmon et al. 
(1987) . Tzeng and Chen (1988) showed that the divergent center over 
the western Pacific Ocean at 200 mb shifted northeastward in response 
to the removal of the mountains, while the northward shift of 
precipitation variance associated with the divergent center is also 
evident in Figure 28b. 

Blackmon et al. (1987) noted substantial differences between the 
precipitation simulated with and without mountains in several areas of 
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the tropics; e.g. , east of the Andes, over the Atlantic and Indian 
Oceans and, most notably, near the East Indies. This is also true in 
the present study when we compare the precipitation rate of the 
control run (Figure 6a) with the no-mountain run (Figure 28a) . If one 
views this as an anomaly experiment, the mountain minus no-mountain 
precipitation difference gives marked heating anomalies throughout the 
tropics. Therefore, the no-mountain climates can not be interpreted 
as sinply due to the lack of orographic effect on the atmospheric 
flow, since the insertion of the mountains has seme effect on the 
distribution of diabatic heating as well. 

Figures 29a and 29b show latitude-time sections of z with m=l and 
eastward periods of 15-20 days at 200 mb and precipitation rate with 
m=l-15 and the period of the Kelvin wave. The two episodes of the 
Kelvin wave are evident in Figure 29a. It is seen that the 
latitudinal width of the Kelvin wave agrees well with that calculated 
from the linear wave theory. Hay ash i and Golder (1978) conducted 
several experiments with a GFDL QCM from which topography, midlatitude 
disturbances and condensational heating were eliminated one by one. 
Based upon space-time spectral analysis of the simulations, they found 
that Kelvin waves are hardly affected by the topography and the 
midlatitude disturbances. 

Similar to that in the control experiment, the band-filtered 
precipitation rate (Figure 29b) oscillates nearly in phase with the 
geopotential height. Comparing these two figures, we find that well- 
defined Kelvin wave is accompanied by the occurrence of the band- 
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Figure 29. Same as Figure 7 , except for the eastward moving periods of 15-20 
days in the no-mountain experiment. Contour intervals is 3 m in 
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filtered precipitation peak near the equator. Therefore we conclude 
that removing mountains from the model has little effect on the 
generation of the Kelvin wave. 

4. Maintenance of tiie Kelvin wave 

Figure 30 presents the energetic budget analysis for EMWs with 
m=l at 200 mb. The kinetic energy of the EMWs is shown in Figure 30a. 
Clearly, there is a concentration of the kinetic energy in the period 
range of 10-60 days, while the maximum kinetic energy with a period of 
20 days has a symmetric structure about the equator and decays away 
from the equator. 

The conversion of available potential energy into kinetic energy 
(Figure 30b) for the QWs with a period of about 12-30 days is 
significant, while appreciable positive conversion taking place 
between 5S to 7N indicates that the conversion of available potential 
energy generated by the latent heat release is an energy source for 
the maintenance of the Kelvin wave. 

Similar to the conversion term, the vertical convergence of 
energy flux (Figure 30c) for the EMtfs with a period of 12-30 days is 
significant between 5S and 8N. 

Figure 30d shows the meridional convergence of wave energy flux. 
It is evident that the wave energy flux associated with the EMW in the 
period range of 12-30 days is transported poleward near the equator 
(diverges) and converges somewhere near 10° latitude aw ay from the 
equator. This feature is consistent with that of the Kelvin wave in 
the control experiment. 


109 


Figure 30. Same as Figure 9, except for the no-mountain 

experiment. Contour intervals are 2 J in (a) ; and 10” 5 
W. in (b), (c), (d) and (e) 
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Hie nonlinear kinetic energy transfer is very weak over the whole 
spectrum except for the period band 15-20 days, with which the Kelvin 
wave gains kinetic energy between 15S and 15N through nonlinear 
interaction. 

From the previous discussion, we find that eliminating mountains 
from the model hardly changes the maintenance of the Kelvin wave, 
except for the conversion of available potential energy into kinetic 
energy, which becomes important in the wave maintenance due to the 
elimination of mountains. 

C. The No-cloud Experiment 

The above analysis suggests that latent heat release associated 
with cumulus convection may be responsible for the generation of the 
Kelvin wave, since this wave appears even in the absence of mountains. 
In order to confirm this possibility, cumulus convection was 
eliminated from the model. As a result, a weakened Kelvin wave and a 
free Rossby wave were identified in the winter season. 

1_. Winter zonal-mean wind 

Convective storm systems can transport meteorological properties 
vertically between the lower and upper troposphere and are important 
in maintaining the large-scale balance of angular momentum, energy, 
and moisture in midlatitudes. Nincmiya (1971) suggested that release 
of latent heat in convection systems is the primary factor that 
intensifies the upper- level large-scale flow through strengthening the 
horizontal temperature gradient. Since cumulus convective does not 
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exist in the no-cloud experiment, the vertical exchange of energy and 
moisture is significantly reduced. Figure 31 presents a height- 
latitude section of zonal mean wind averaged over the winter season. 

As expected, due to the restraint of latent heat release, radiative 
effects, surface and other energy balances, the intensity of the jet 
streams is decreased markedly and the jet cores are located more than 
2 km lower than the height where they were situated in the control 
experiment. 

2. Spectral analysis 

Figure 32a shows power spectrum of u at 200 mb and 2S as a 
function of wavenumber and frequency. It is seen that P(u) is more 
spread out in the wavenumber-frequency domain than in the control 
experiment. There are two spectral peaks at m=l and eastward periods 
of 10-15 days and westward periods of 17-24 days. 

Figures 31b and 31c display P(z) and P(u) associated with m=l at 
200 mb. It is evident that the power spectra EMWs are dominated by 
disturbances which are symmetric with respect to the equator and decay 
away from the equator. On the other hand, P(z) for the WMWs is 
dominated by disturbances from extratropical areas in both 
hemispheres, while P(u) for the WMWs has an appreciable peak at 
periods of 20-30 days on the equator. We also notice that the 
latitude-frequency power spectrum of v with m=l does not have 
appreciable spectral peaks in the frequency bands corresponding to 
where spectral peaks of P(u) are observed (not shown) . 

Shown in Figures 31d and 31e are the vertical structures of P(z) 
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over the northern winter. Contour intervals are 1.5 
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and P(u) with m=l at 2S, respectively. It is evident that P(z) is 
essentially dominated by EMWs. There are two spectral peaks with 
eastward periods of 6 days and 20 days locate above the 200 mb level. 
No significant spectral peak occurs at the period range of 10-15 days. 
However, P(u) is dominated by the two spectral peaks which are 
observed in Figure 31c at 200 mb. 

3. Propagation characteristics 

From the dispersion diagram associated with m=l of the Laplace's 
tidal equation (Longuet-Higgins, 1968) , we can estimate the equivalent 
depth of the spectral peak associated with the WMW with m=l and a 
period of 17-24 days. For this peak, there are two possible locations 
on the dispersion diagram. Hie location may be either on the n=l or 
the n=3 curve, where n is the latitudinal nodal number. Thus the wave 
associated with this peak can be interpreted as either 1) an 
equatorially trapped Rossby wave with n=l, or 2) a free Rossby wave 
with n=3. If the first interpretation is correct, the meridional 
velocity associated with the wave should reach a maximum at seme 
latitude, which is estimated to be 17°, according to the dispersion 
diagram. However, the latitudinal distribution of P(v) (not shown) 
increases poleward without showing a clear indication of the v 
maximum. 

If we adopt the second interpretation, the estimated equivalent 
depth (8.4 km) is very close to the equivalent depth of free 
oscillations, and the gradual poleward decrease of P(u) and the 
monotonic increase of P(v) from the equator to 17° latitude are more 
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consistent with the theoretically predicted structure of the wave 
(Longuet-Higgins, 1968) . The existence of the second gravest 
symmetric free Rossby mode with a period of 16-18 days has been 
reported by many authors [see Madden (1979) and Walkerscheid (1980) 
for reviews]. Therefore, this peak is identified as the second 
gravest symmetric free Rossby mode with n=3 . In the following 
discussion we shall focus on the spectral peak associated with m=l and 
periods of 10-15 days. 

The spectral analysis indicates that the disturbance associated 
with m=l and eastward moving periods of 10-15 days is a Kelvin wave. 
Figure 33 shows the horizontal distributions of geopotential height 
and wind components of the disturbance at 200 mb. The most striking 
feature is the reduction of both the wave amplitude and latitudinal 
extent, compared to those in the control experiment. The amplitude of 
geopotential height and zonal wind are about 3 m and 3 ms - -*-, 
respectively. The latitudinal width is about 5° latitude. However, 
based upon linear wave theory which assumes the atmospheric Kelvin 
wave is a free wave, the latitudinal width and the vertical wavelength 
are estimated to be 16° latitude and 11 km (Table 3) , respectively. 

The discrepancy between the model simulation and the linear theory 
predication for latitudinal width strongly suggests that the 
atmospheric Kelvin waves are forced waves instead of free wave modes. 

Figure 34 shows a latitude-time section and a height-time section 
of the geopotential height oscillations with m=l and periods of 10-15 
days, repectively. In Figure 34a, the Kelvin wave with reduced 
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KELVIN WAVE (m=l) FEB 21 



Figure 33. Same as Figure 4, except for the Kelvin wave with m=l 
and periods of 10-15 days in no-cloud experiment. 
Contour interval is 1 m 
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latitudinal extent is evident during days 70-110 and after day 135. 

Our results confirm that obtained by Hayashi and Golder (1978) who 
showed that the power of the Kelvin wave diminished drastically after 
condensational heat as well as midlatitude disturbances, was 
eliminated from the GCM they used. Thus, it is revealed that Kelvin 
waves may be generated due to the effect of latent heat release 
produced by cumulus convection in the troposphere. 

The vertical structure of z at 150E longitude and 2S (Figure 34b) 
reveals a markedly downward phase progression accompanied by the 
appearance of the Kelvin wave (days 70-120) . A noticeably reduced 
amplitude is clear, implying weaker upward energy propagation. Since 
clouds are not formed in this experiment, the solar radiation can 
reach the earth surface without being reflected back to space by the 
clouds. Therefore, the atmospheric circulations will be modified by 
the incoming solar radiation. 'Rseng and Chen (1988) noted that three 
pronounced divergent centers in the no-cloud experiment appear at 
about 10S over the three tropical continents (the American, African, 
and maritime continents) in response to the landmass heating. 

However, thermally direct convections associated with the three 
divergent centers are not able to transport sufficient energy upward 
to stimulate the Kelvin wave in the absence of the cumulus convective 
process. 

£. Maintenance of the Kelvin wave 

Figure 35 presents energetic budget analysis for EMtfs with m=l at 
200 mb. It is evident that the kinetic energy (Figure 35a) for those 
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EMWs with periods of 10-30 days is significant. The maximum kinetic 
energy with a period of 10-15 days appears on the equator with an even 
structure about the equator. 

The conversion of available potential energy to kinetic energy 
shown in Figure 35b is significant for the 0W with a period of about 
6 days. However, over the period band of 10-15 days, with which the 
Kelvin wave is observed, the conversion is of less importance for the 
maintenance of the Kelvin wave. Only a weak positive conversion lying 
between the equator and 5S indicates less energy input into the Kelvin 
wave than that in the control experiment. 

As expected, eliminating cumulus convection from the model will 
drastically suppress vertical transport of energy and moisture. It is 
not surprising to see that the vertical convergence of wave energy 
flux (Figure 35c) is almost negligible due to the reduction of 
vertical energy flux. Therefore, it reveals that the two primary 
energy sources of the Kelvin wave in the control experiment, the 
conversion of available potential energy to kinetic energy and the 
vertical convergence of wave energy flux generated by the cumulus 
convections in the tropics, are no longer able to provide sufficient 
energy to the Kelvin wave in this case. 

In contrast to that in the control experiment, the Kelvin wave in 
the no-cloud case is not accompanied by meridional divergence of the 
wave energy flux due to the pressure work (Figure 35d) . Weak 
divergence and convergence with a period of 12-15 days located on 
either side of the equator indicate that there is a weakly southward 


Figure 35. Same as Figure 9 , except for the no-cloud experiment 

over the northern winter. Contour intervals are 1 J in 
(a); and 10“ 5 W in (b) , (c) , (d) and (e) 
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Figure 36. Same as Figure 10, except for the no-cloud experiment 
over the northern winter. Contour interval is 0.5x10' 
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energy flux near the equator associated with the Kelvin wave. 

As shown in Figure 35e, the nonlinear energy transfer is 
significant for the EMWs with periods of 4-6 days and 8-30 days. 
Separating the nonlinear energy transfer into wave-mean flow and wave- 
wave interaction, we find that the suppressed zonal flow weakens the 
wave-mean flew interaction as shown in Figure 36a. It is noteworthy 
that there is a weak kinetic energy input into EMWs with periods of 
10-30 days north of 5N. On the other hand, enhanced energy transfer 
due to the wave-wave interaction is evident over the whole spectrum 
(Figure 36b) . It is noted that EMWs with periods of 8-15 days, in 
which the Kelvin wave is observed, gain energy north of 5N and between 
the equator and 10S. 
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VI. OONCLCDING REMARKS 

The equatorial planetary waves simulated by the NCAR CCM1 general 
circulation model have been investigated in terms of space- time 
spectral analysis (Kao, 1968; Hayashi, 1971, 1973) and energetic 
analysis (Hayashi, 1980) . These analyses are particularly applied to 
grid-point data on latitude circles. In order to test some physical 
factors which may affect the generation of tropical transient 
planetary waves, three different model simulations with the CCM1 (the 
control, the no-mountain and the no-cloud experiments) have been 
analyzed. 


A. Summary and Discussion 

On the basis of the analyses of the three experiments, the 
identified equatorial planetary waves can be summarized as following: 

1. two Kelvin waves are identified in the summer of the control 
experiment, one at m=l and an eastward propagating period of 
18-25 days with a vertical wavelength of about 5.9 km, and 
the other at m=2 with an eastward moving period of 5-6 days 
and a vertical wavelength of 12 km. On the other hand, in 
the winter of the control experiment, the Kelvin wave with 
m=2 is absent and only the Kelvin wave with m=l and a period 
of 15-20 days is found. The characteristics of these Kelvin 
waves are in agreement with observations by Zangvil and Yanai 
(1980) and Yanai and Lu (1983) . 

2. Kelvin waves at m=l and with a period of 15-20 days, and at 
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m=l and a period of 18-25 days are clearly identified in the 
summer and winter of the no-mountain experiment, 
respectively. The amplitudes and propagation characteristics 
of the Kelvin waves are hardly affected by the removal of 
mountains, compared to those in the control experiment. 

3. lhe Kelvin wave at m=l and with a period of 10-15 days is 
detected only in the winter of the no-cloud experiment. 
However, the amplitude and the north-south extent of the wave 
are significantly reduced, in agreement with the results 
obtained by Hayashi and Golder (1978) . 

4. A mixed Rossby-gravity wave is found with periods of 12-15 
days for m=5 in the summer of the control experiment and has 
a vertical wavelength of about 3.3 km. flie amplitude of the 
meridional wind associated with the wave is about 1 ms -1 . 
However, observations (Zangvil and Yanai, 1980; Yanai and 
Lu, 1983) reported the existence of a mixed Rossby-gravity 
wave with m=4 and a period of 5 days at 200 mb. On the other 
hand, this wave disappears completely in both the no-mountain 
and the no-cloud experiments. 

5. A equatorially trapped Rossby wave associated with a 20-30 
day period for m=3 and a vertical wavelength of about 14 km 
is observed in the winter of the control experiment. 

However, observed equatorially trapped Rossby waves are found 
at m=2 and with periods of 10 days and 30 days, respectively. 
'This wave is absent in both the no-mountain and the no-cloud 
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experiments. 

6. In addition, the next gravest symmetric free Rossby wave at 
m=l and with westward period of 17-24 days is detected in the 
no-cloud experiment. 

Based upon the energetic analyses, the physical processes which 
are important in maintaining the identified equatorial waves in the 
present study are showed in Table 4. Hie maintenance of the waves can 
be sutimarized as following: 

1. The two primary energy sources for maintenance of the Kelvin 
waves, the mixed Rossby-gravity wave and the equa tori ally 
trapped Rossby wave wave modes at 200 rob are the conversion 
of available potential energy to kinetic energy and the 
vertical convergence of wave energy flux. The Kelvin waves 
are accompanied by a meridional divergence of wave energy 
flux near the equator, while the mixed Rossby-gravity wave 
and the equatorially trapped Rossby wave are associated with 
meridional convergence of wave energy flux. The latter 
finding suggests that these westward moving waves may be 
excited by the pressure work due to midlatitude disturbances, 
i.e. , the lateral forcing. 

2. The two primary energy sources in both the control and the 
no-mountain experiments are weakened markedly in the no-cloud 
experiment. The vertical energy flux is especially 
negligible, while the meridional wave energy flux associated 
with the Kelvin wave is weakly southward instead of converging 



Table 4. Various physical processes associated with wave maintenance 
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on the equator. 

3. Hie nonlinear energy transfer associated with the Kelvin 
waves is of less importance compared to that associated with 
the mixed Rossby-gravity wave and the Rossby wave. Hie 
kinetic energy transfer through the nonlinear interaction 
appears to be more important in the winter than in the 
simmer. 

4. CXir results indicate that the linear waves theory is still 
valid even though the nonlinear effect may slightly change 
the frequencies arid the structures of the equatorial waves. 

Tsay (1974) showed that the Kelvin wave in the NCAR GCM is 
accompanied by meridional energy convergence between 12 km and 18 km 
over the equator. However, our results show that the Kelvin waves in 
the COU are associated with meridional energy divergence near 200 mb, 
which is in agreement with observation by Yanai and Lu (1983) . 

Hie spectral analysis of the model precipitaion rate gives a 
stochastic distribution in the wavenumber-frequency domain. On the 
other hand, there are high correlations between Kelvin waves and the 
band-filtered precipitation rate that is symmetric about the equator, 
and between the mixed Rossby-gravity wave and the band-filtered 
precipitation rate that exhibits an antisymmetric structure with 
respect to the equator. This finding confirms Holton's random thermal 
forcing theory, which states that the frequencies of the Kelvin wave 
and the mixed Rossby-gravity wave are determined by a certain 
filtering effect of the atmosphere and not by the wavenumber-frequency 
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distribution of the heating function. In the same sense , the 
fluctuations of the Kelvin waves and the mixed Rossby-gravity waves in 
the present study are more accurately characterized as episodic, rather 
than steady oscillations (Salty et al. , 1984) . This can be 
interpreted as both waves originating in the stochastic latent heat 
release of troposphere cloud ensembles. 

The band-filtered precipitation rate is highly correlated with the 
Kelvin waves and the mixed Rossby-gravity wave at 150E longitude, at 
which the strongest divergent center is located (Tzeng and Chen, 

1988) , indicating that the cumulus convective process in this region 
is resposible for transporting energy from lower levels to upper 
levels in order to maintain the waves. 

B. Deficiencies of the Model Simulations 

In view of these model simulations, we find that the NCAR CCM1 
can simulate equatorial planetary waves well to a certain extent; in 
particular, the Kelvin waves have been simulated successfully. We 
believe the reasons for this success are rooted the treatment of 
cumulus convection and the different parameterization schemes. 

The difference between the model simulations and observations is 
that the simulated WMWs tends to have smaller horizontal wavelengths. 

A possible factor causing the defect could be the coarse resolution of 
the model (im=15) . 

Canparing with the vertical propagation behavior of Kelvin waves 
above the tropopause (Salby et al., 1984; Hayashi et al., 1984), we 
note that the vertical structures of the Kelvin waves exhibit little 
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vertical phase tilt in the troposphere, but a vertical phase change of 
180 near 500 mb. This discrepancy can be explained in terms of the 
near-field behavior of Kelvin waves as discussed by Salby and Garcia 
(1987) . Based upon a set of linearized primitive equations in 
spherical geometry, they showed that the Kelvin response arising from 
a heating centered on the equator in the upper troposphere is 
characterized by a 180° phase change in the vertical before radiating 
entirely out of the troposphere. 

On the other hand, the model's top boundary may affect Kelvin 
waves significantly. Comparing two model simulations by the CCM with 
different top levels (at 0.1 mb and 10 mb), Boville and Cheng (1988) 
found that Kelvin waves in the model with the top level at 10 mb are 
reflected by the top boundary and become vertically standing 
oscillations which retain approximately their original zonal phase 
speed. 


C. Suggestions for Future Study 
In the present study, we treated equatorial transient waves 
as the fundamental solutions of a zonally symmetric basic state 
predicted by linear wave theory. However, the consideration of 
waves in a nonzero, zonally asymmetric basic state has been proven to be 
important to climate dynamics. Webster and Chang (1988) proposed that 
the longitudinal asymmetry of the basic flew in the tropics so alters 
the characteristics of equatorially trapped modes that their 
Doppler-shifted group speed becomes locally zero, producing a local 
convergence of energy to the east of the equatorial westerlies. We 
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plan to extend our study by taking into account zonally asymmetric 
states. By doing this, we may be able to investigate the wave 
transition between the region of easterlies and westerlies along the 
equator, and test whether or not the transition occurs as an evolution 
of the latitudinal modal structure of the equatorial modes as the 
modes propagate from an obstensibly easterly regime to a energy 
accumulation point within the westerlies. 
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Figure A2. Same as Figure 6 , except for the northern win 1 
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Figure A3. Same as Figure 7, except for the eastward moving 
periods of 15-20 days over the northern winter 
































